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SUMMARY 


Great progress has been made in recent years toward improving 
the reliability, speed, and accuracy of the various experimental 
test procedures available to the aircraft designer. This paper 
describes the development of such a procedure whereby flying 
scale models, which are dynamically similar scale duplicates of 
the full-scale prototypes, may be accurately flight-tested through 
the use of positioning, multichannel radio remote control. This 
development allows the designer to accurately determine many 
dynamic characteristics of a new design which would not or- 
dinarily be available until the actual airplane is flight-tested. 
This procedure is particularly well adapted to the determination 
of dynamic stability during the take-off and landing of flying 
boats and was developed primarily for this purpose. However, 
it is directly applicable to many other dynamic problems of 
both landplanes and seaplanes. 

The design and construction of the self-propelled free flight 
models and the remote control equipment are considered, as 
well as a general discussion of the instrumentation, analysis, and 
testing procedure. It is concluded that: 

(1) Full-scale hydrodynamic stability can be accurately 
predicted by scale, free body, dynamically similar models through 
the use of multichannel radio remote control. 

(2) The special Consolidated Vultee radio transmitter and 
receivers, designed for the proportional control of seven simul- 
taneous channels, have proved entirely feasible and reliable. 

(3) The special two-cylinder, two-cycle gasoline engines have 
excellent speed control characteristics and have proved satis- 
factory and dependable under the most severe conditions of salt 
water and continuous high-power operation. : 

(4) The method of instrumentation and analysis for free 
bodies has proved to be simple, rapid, accurate, and extremely 
Versatile in its operation. 

These conclusions are substantiated by extensive operation 
and the results of full-scale correlation. 


, 
INTRODUCTION 


NE OF THE most difficult problems in aircraft de- 
sign is the accurate prediction of the complex 


Variables involved in the dynamics of free bodies in 
Motion. Perhaps the best example of this problem is 
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the stability phenomena associated with the take-off 
and landing of a flying boat. The most satisfactory 
approach to this problem has been the development and 
use of dynamically similar models? which are not only 
geometric scale duplicates of a contemplated or exist- 
ing full-scale airplane but are also to scale in gross 
weight, moment of inertia, power, and aerodynamic 
forces. This resolves the analysis into a direct experi- 
mental procedure inasmuch as the model can be flown, 
and being dynamically similar, it will reproduce ac- 
curately the full-scale motions as well as the rates of 
movement. 


The handicap to this form of analysis in the past has 
been the difficulty involved in the control and instru- 
mentation of the free bodies, which invariably led to 
partial restraint of some nature. Attempts have been 
made in the towing basin to give dynamic models of 
seaplanes freedom about all three axes, and even fore- 
and-aft freedom, but the results have not proved fully 
satisfactory when any more than freedom in rise and 
pitch is allowed. Even here the rise is so limited that 
the models can no more than just take off and the fore- 
and-aft acceleration is that impressed by the rigidly 
connected towing carriage. In addition, all forms of 
partial restraint restrict all sideways motion, preclud- 
ing any investigation of dynamic directional stability. 
For dynamic studies of landplanes the problem is even 
more difficult inasmuch as there is no equipment avail- 
able similar to the seaplane towing carriage for towing 
the models and the free flight wind tunnels are seri- 
ously handicapped by size and maneuverability. 


To offset these disadvantages and to obtain unre- 
strained free body operation, the Consolidated Vultee 
Hydrodynamics Department, in cooperation with 
members of the Radio Laboratory, started a project 
late in 1943 to develop free body, dynamically similar 
flying models to be remotely controlled by positioning, 
multichannel radio. 
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Three views of the radio model. Top, front view; 
center, plan view; bottom, side view. 


Fic. 1. 


In order that the maximum benefit would result from 
the experimental project, it was decided that the first 
free model should be one that had ample towing basin, 
wind-tunnel, and full-scale test results to which the 
free flight model test results could be compared. Also, 
the first experimental model should be simple, reducing 
as far as possible the complexity of the problem. With 
these considerations in mind, the Consolidated Vultee 
Model 31, a twin-engined flying boat, appeared to be 
the logical choice. To make the project possible it 
was necessary to develop an entirely new principle of 
radio transmission and reception; small lightweight 
positioning servomotors; lightweight, dependable gaso- 








Fic. 2. The full-scale Model 31 airplane. 
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line motors; and a simple but accurate means of re. 
cording and analyzing free body test data. 


MopDEL DESIGN 


The model used in this development (Fig. 1) is a 
one-eighth scale, dynamically similar model of the 
Consolidated Vultee Model 31 twin-engined flying 
boat shown in Fig. 2. It is similar to the N.A.C.A. 
model reported in reference 1 and shown in Fig. 3a, 
The main points of difference are: 

(1) The N.A.C.A. model was powered with two 
Sawyer variable frequency, alternating-current, induc- 
tion motors rated at 4,700 r.p.m. for scale thrust; 
whereas the radio model is equipped with two, two- 
cylinder, two-cycle gasoline motors rated at 4,250 
r.p.m. 

(2) The N.A.C.A. model was equipped with a twin 
tail similar to the tail assembly on the full-scale air- 
plane; whereas the radio model is equipped with the 
new single tail that had been proposed for the full- 
scale production airplane, the wind-tunnel model of 
which is shown in Fig. 3b. 





Fic. 3b (bot- 
The Production Model 31—wind-tunnel model. 


Fic. 3a (top). The N.A.C.A. towing basin model. 


tom). 


(3) To increase the lateral stability and facilitate 
control of the free radio model, a slight amount of 
additional dihedral was added to each outer panel. 
As the N.A.C.A. model was restrained in yaw and roll 
this was not required. 

The model is constructed of '/,9-in. skin, balsa mono- 
coque planked onto '/s-in. plywood bulkheads, as 
shown in Fig. 4. The entire exterior surface is stripped 
with rice tissue paper laid ir shellac which seals the 
porous balsa and provides a glossy finish. Three coats 
of spar varnish on the tissue provide a light, watertight 
construction. 

The wing is equipped with full-span, dural leading- 
edge slots which compensate for scale effect on lift due 
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to Reynold’s Number. Aerodynamic tests were con- 
ducted on the model prior to the stability tests (refer- 
ence 3), and it was established that the slots resulted in 
full-scale slope to the lift curve and provided full-scale 
Clmaz- ‘Lhe equipment for conducting these tests is 
shown in Fig. 5 where six component data were ob- 
tained through the use of an electric strain gage bal- 
ance. Comparison with wind-tunnel results showed 
good agreement of lift, moments, and control effective- 
ness. The wing is also equipped with scale Fowler 
faps that are remotely controlled between 20° and 
49°, Provision for 0° flaps was not made because of 
fap track complication. All control surfaces are 
covered with doped silk and are remotely controlled by 
radio. 

The model has a wing area of 16.38 sq. ft. and a span 
of 13.75 ft. The weight of the model complete with 
engines, propellers, fuel, radio receiver, servos, and 
batteries is 78 lbs., or 40,000 Ibs. full scale. Since the 
scale gross weight is approximately 100 Ibs., this per- 
mits 22 lbs. of lead ballast to be added to balance the 
model with the pitching moment of inertia to scale. 





Radio model construction view. 


Fic. 4. 


POWER PLANT 


The development of a self-contained power plant 
that would have sensitive speed control, light weight, 
and rugged dependability posed a particularly difficult 
problem. In the case of free bodies it was not possible 
to use the special variable frequency electric motors 
that had been developed for the partially restrained 
towing basin models. After careful study it was de- 
cided to develop a special two-cycle gasoline engine 
for the job. Power-plant specifications were prepared 
and Ohlsson and Rice Manufacturing Company, of 
Los Angeles, agreed to make the engines, which were to 
fit into a 9-in. cowl, weigh approximately 3 Ibs., de- 
liver 11/2 b.hp. at 5,000 r.p.m., and have complete 
throttle control from idling to full throttle. 





Fic. 5. The aerodynamic electric strain gage test equipment, 

Four months after the contract was signed the first 
experimental engine was completed and bench tested. 
The engine is two-cylinder two-cycle and has a weight 
of 3.5 Ibs. It is nominally rated at 1.6 b.hp. at 4,250 
r.p.m. which gives scale power and r.p.m. for the R- 
3350 engine on the full-scale airplane. Dynamometer 
tests using a one-eighth scale duplicate of the full-scale, 
three-bladed, 16-ft. Curtiss Electric propeller, resulted 
in producing scale static thrust at scale r.p.m. 

The power-plant installation is shown in Fig. 6 and 
weighs 8.0 Ibs. complete per nacelle. This weight in- 
cludes the engine, propeller, coil, gasoline tank (which 
is located in the top of the nacelle just behind the fire 
wall), engine mount, fire wall, servomotor, and cowling. 
The engine has a displacement of 2.77 cu. in. and an 
overall span of 9.5 in. including the spark plugs. The 
cylinders fire simultaneously and the crankcase is 
supercharged by a direct drive rotary impeller. The 
induction system includes a manual mixture control 
for starting and an automatic choke, butterfly, and 
spark control for varying the mixture while running 
at varying speeds. This feature provides exceptional 
throttle control for a two-cycle engine. The intake 





Fic. 6. Radio model power-plant installation 
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valve is rotary and the exhaust is valved by the piston. 
The propeller shaft and hub are splined with standard 
S.A.E. splines. The propellers are machined from 
24ST dural and are adjustable pitch. 

During the early stages of development, considerable 
trouble with engine stoppage was experienced. This 
was caused primarily by faulty carburation and salt 
spray shorting the ignition system. During the last 80 
runs this has been corrected by using direct fuel injec- 
tion and redesigned shielding, which has resulted in 
only four engine failures. In most cases this is only 
one engine and control is maintained. The model has 
a gasoline capacity for 25 min. of operation and, ex- 
cept for the cases of engine or radio failure, the model 
is beached under its own power at the end of each run. 
This is accomplished by taxiing to the ramp and the 
model is slowed down by intermittently cutting the 
ignition, which is remotely controlled from the trans- 
mitter. When a few feet from the ramp, the ignition is 
out completely and the model is floated into its cradle 
and pulled up onto the ramp for refueling or changing 
of the center of gravity. 


SERVO SYSTEM 


The elevators, rudder, ailerons, flaps, and independ- 
ent throttles are operated by small, electric positioning 
servomotors as illustrated in Fig. 7. These servos are 
five-pole, split field electric motors geared to a jack- 
screw and rated 1 amp. at 6 volts. The stalled thrust 
on the jack is 14 lbs. All servos are equipped with 
potential dividers that position the servo proportion- 
ately to the modulated audio signal. Each servo con- 
tains a homing circuit and, in case of emergency, the 
transmitter is turned off and all controls on the model 
home to a predetermined position of gliding flight with 
idling propellers. If the pilot so desires, he may switch 
the transmitter on again and pick up control of the 
model. This is also a safety feature that automatically 


Sidi Jd 


FULL SCALE 





Fic. 7. Radio model electric servomotor. 





The seven-channel radio transmitter. 


Fic. 8. 


comes into action should the transmitter fail for any 
reason. To prevent overriding and hunting the ser- 
vos are equipped with magnetic armature brakes that 
are instantly released whenever an operating current 
flows to the servo. 

For starting and ground checking the equipment, a 
booster system has been designed. The booster box 
contains storage batteries and leads connected to a 
switching panel. The model contains a boost plug-in 
so that the ship’s batteries may be cut out and all 
power is drawn from the booster batteries. Each 
servo is equipped with jack plugs so that the leads from 
the booster switching panel may be plugged into any 
servo and that servo can be checked or operated inde- 
pendently of the radio. For starting, the booster 
leads are plugged into the throttle servos and the 
throttles are manually controlled until the engines are 
adjusted and the warm-up is completed. The booster 
leads are then pulled and the power lead is disconnected, 
which automatically cuts in the ship’s batteries. The 
radio switch is turned on and the pilot at the trans- 
mitter takes over and checks out all controls on the 
beach. At the pilot’s signal the model is launched into 
the water and the test run is started. 


RADIO CONTROL 


The radio equipment was designed and built by the 
Consolidated Vultee Radio Laboratory in close co- 
operation with the Hydrodynamics Department. The 
transmitter, shown in Fig. 8, is built in the form of a 
full-scale flight deck with adjustable pilot’s, seat and 
full-scale flight controls, and is highly mobile. The 
transmitter is four-stage, crystal-controlled, and has a 
63-watt output. The single R.F. carrier carries seven 
audio frequencies corresponding to the seven primary 
remote-control channels. This allows simultaneous 
and independent control of flaps, ailerons, rudder, 
elevator, two throttles, and ignition. The position of 
any control is determined by the modulated amplitude 
of the corresponding audio frequency. 
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Fic. 9. The seven-channel radio receiver. 


The transmitter is equipped with trimming knobs 
that trim the zero positions of all the channels. This 
feature saves a considerable amount of time that might 
otherwise be spent in adjusting for the slight day-to- 
day drift in the synchronizing of the radio equipment. 
Inasmuch as the positioning feature of the control 
system is predicated on amplitude modulation of the 
audio signal, an extremely sensitive constant volume 
control was necessary. The transmitter will maintain 
constant signal strength up to a range of 3,500 ft., 
which is beyond the practical range of sight for satis- 
factory ground control. If the model exceeds this 
range the maximum amplitude of the signal will grad- 
ually reduce, which is coordinated in such a manner 
that the model goes into a gradual left turn and returns 
into range. 

The pilot’s control wheel is equipped with a micro 
thumb switch that engages the autopilot. It is possible 
to fly the model manually and, at any time, at the press 
of this switch, to engage the autopilot to hold any 
stabilized position. With the autopilot engaged the 
rudder pedals are disengaged and the pilot’s wheel 
controls the model through the gyro. A rate gyro in 
the model operates the rudder automatically to produce 
coordinated turns as determined by the aileron control. 
For the making of high-speed stability runs on the water 
where independent rudder control is desired, a switch 
on the instrument panel of the transmitter may be 
thrown, which leaves the rudder full manual when the 
autopilot is engaged. Through the use of this circuit 
it is possible to stabilize pitch and roll while maintain- 
ing full manual rudder control. 

The radio receiver is shown in Fig. 9 and consists of 
the pickup, filtering, modulation, and relay circuits. 
This receiver, including the waterproof dural case, 
weighs 12.5 Ibs. complete. The changes in modulation 
of the audio for any channel operate a sensitive micro- 
relay with the plate voltage. This relay, in turn, 
operates a power relay that carries the servo current. 
The positioning potentiometer on the servo balances 
the bridge containing the plate voltage and provides the 
control follow-up and positioning feature. Through 
the use of magnetic braking on the servo armature, 
which eliminates overriding and hunting, the sensitivity 
is of the order of one part in 40, which amounts to less 
than 1° on most of the controls. 


Fig. 10 shows the complete receiver and battery in- 
stallation in the hull of the model. In this figure the 
receiver is shown in its waterproof dural box. The 
total weight of the receiver and batteries is 24.5 Ibs. 
Both units are clamped to longitudinal dural tracks 
and may be moved fore and aft to secure any c.g. posi- 
tion within the range of 10 to 40 per cent of the mean 
aerodynamic chord. All electrical connections are 
terminated in quick disconnect gang jack plugs so 
that the receiver may be readily removed and no errors 
can be made in reconnecting the elaborate circuits. 
On the top of the battery box are located the two igni- 
tion cutoff relays. These relays work direct without a 
servo. The battery box contains ‘‘B” and ‘‘C”’ bat- 
teries for the radio, plus two small 6-volt and two 2- 
volt storage batteries for power. These storage bat- 
teries have plastic cases and are nonspillable. The 
batteries are arranged into two 8-volt circuits totaling 
14 amp.-hours’ output. One circuit is for the ignition, 





Fic. 10. Radio receiver and battery installation. 





WATER SPEED INDICATOR 


Fic. 11. Cutaway inboard profile of the radio model. 
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which is kept separate from the radio and servo func- 
tions. Fig. 11 shows a cutaway inboard profile of the 
model from which the relative locations of the various 
pieces of equipment may be obtained. 

The dependability of the transmitter, receiver, and 
servos is excellent. It is customary to check the radio 
equipment every 30 runs, which is calculated to repre- 
sent 30,000 independent operations. During the train- 
ing, familiarization, experimentation, and testing of 
this model, over 150 runs were completed, and there 
were only three radio failures. These were usually 
minor in nature and consisted of losing one or two 
channels because of vibration or occasional salt-water 
spray. During the last 80 runs there have been no 
failures. 


INSTRUMENTATION AND ANALYSIS 


One of the greatest problems encountered in de- 
veloping an accurate, rapid method of testing free 
bodies by remote control was the development of a 
sound but simple system of instrumentation and anal- 
ysis. Original efforts were directed toward equip- 
ping the model with test instruments that were recorded 
photographically in the model. Toward this end, a 
small 8-mm. motion-picture camera was developed 
which weighed but 5 lbs. and was mounted on the flight 
deck. This camera recorded air speed, trim, time, and 
a manometer reading that indicated the moment of 
leaving the water. An N.A.C.A.-type of visual trim 
indicator was built into the windshield of the model, 
which projected an image of the horizon onto a ground- 
glass grid which was photographed. In conjunction 
with the trim indicator, a small, sensitive, pitot-ven- 
turi air-speed indicator was developed and mounted in 
the field of the camera. A stop watch and manometer 
connected to the hull bottom at the step completed 
the installation. 

This method was checked extensively but indicated 
several severe shortcomings. The camera was con- 
nected to the throttle and would last only one run 
without rewinding, necessitating returning to the 
beach, and the air-speed indicator did not produce the 
accuracy considered necessary, particularly in the 
region of the hump. To offset these disadvantages, 
there was revived a method that had been developed 
many years ago by Consolidated Vultee to photogra- 
phically analyze full-scale, rough-water tests‘ without 
the use of elaborate instrumentation. This method 
has proved to be extremely simple and accurate, and 
has been adopted as standard procedure. The instru- 
mentation and procedure will be briefly outlined. 

The only instrumentation required for this method 
of analysis is a good 35-mm. motion-picture camera 
with a telephoto lens, mounted on a tripod, and a 
special film analyzer designed by the Hydrodynamics 
Group and shown in Fig. 12. This method immed- 
jately eliminates the element of mechanical failures in 
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Fic. 12. Step-by-step film analysis procedure, using the Hydro- 
dynamic Film Analyzer. 





MOTION PICTURE 
CAMERA 





SEAGOING LAUNCH 











ae 


Fic. 13. Photographic instrument arrangement. 
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Fic. 14. Geometry determining the true angle of trim. 
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elaborate instrumentation, with the accompanying 
delays when tests have to be rerun. 

Since this system is based on a theorem of descrip- 
tive geometry dealing with the relation and propor- 
tionality of true length to foreshortened length, it is 
necessary only to determine a fixed longitudinal and 
lateral reference on the model, the true lengths of which 
are known, and to record on photographic film the 
foreshortened relationship of these lengths from some 
point away from the model. As it is necessary only 
to obtain the proportionality, the distance from the 
camera to the model need not be known and can even 
vary during a run. Therefore, the camera is set up 
beside the transmitter and takes a picture of the model 
as the pilot sees it, always keeping the horizon, which 
is the trim reference, in the picture. 

Two large black crosses are painted on the hull of 
the model as shown in Fig. 22. These crosses are par- 
allel to the deck line and, being a known distance apart, 
constitute the longitudinal reference dimension. As 
the propeller discs or nacelles are visible in all runs, the 
distance between the propeller axes is used as the 
lateral reference dimension. The analysis is made by 
the motion-picture record by frames in the film an- 
alyzer. The camera must run at constant speed so 
that each frame covers a definite and known incre- 
ment of time. 

Because of the relative motion between the camera 
and the model, it is first necessary to correct for angu- 
larity between the camera and the path of the model. 
This is illustrated in Fig. 13. As the frames were run 
through the analyzer, the foreshortened values of A 
and B, which are the longitudinal and lateral reference 
dimensions, respectively, were tabulated. The tan- 
gent of the angle, y, will therefore be, 


tan Y = (A/B) X (True Length B/True Length A) 


from which y can be determined. This fixes the path 
angle of the model to the camera. 

The foreshortened angle of trim, 6, is then read by 
measuring the angle between the crosses and the 
horizon, This is based on the assumption that the 
horizon is infinite and provides a true horizontal. This 
angle can be corrected for the angularity, y, to give the 
true angle of trim, 7. The geometry of this correction 
to the foreshortened angle is illustrated in Fig. 14. 
Faired values of y and 8 are used to determine 7. In 
Fig. 14, 


sin y = L,/L 
L = L,/sin y 
and 
sin 8B = D/L, 
D = L, sin B 
also 
sinr = D/L 


= L,sin 8 sin y/L, 
therefore, 


sin r = sin y sin 8 


From this relationship, trim versus frame or time is 
plotted as in Fig. 15. 

It can be noted in this sample plot that readings were 
taken every two frames. As the camera speed was 24 
frames per sec., this is a reading of trim every 0.08 sec. 
Obviously, such accuracy is not usually necessary and 
it is customary to analyze only every tenth frame, or a 
reading every 0.417 sec. If necessary, and in doubtful 
cases, it is possible to read every frame, which is an ac- 
curate reading every 0.04 sec. As a rule, a run re- 
quires approximately 10 sec.; therefore, a complete 
analysis can be secured by running 24 or 25 frames 
through the analyzer, with possibly two or three extra 
frames at critical points. 

The run shown in Fig. 15 was run on a particularly 
rough day and is presented as a sample to show the 
accuracy with which trim can be determined. It can 
be noted that in the slow speed region up to hump speed, 
it is possible to plot the actual effect of wave action on 
the model as it moves through the rough water. This 
run was carried to a speed where the model was leaving 
the water at the top of the wave crests. 

The speed record is also obtained from every tenth 
frame. Here, too, greater accuracy can be secured if 
deemed necessary. By observing the speed plot in 
Fig. 15, however, it can be noted that every tenth frame 
gives sufficient accuracy to actually see the slowing 
down through the hump. The speed is determined by 
computing the average speed between frames. This is 
done by measuring the progressive distance the image 
of the model moves across a fixed point on the horizon. 
For this purpose buoys are fixed in the background to 
be used as speed references. If from frame 200 to 
frame 210 the front cross on the hull image moves from 
3.50 in. to 6.85 in. ahead of the fixed speed reference 
point, the model has moved 3.35 in. in 0.417 sec. If 
the average longitudinal reference distance (the dis- 
tance between the crosses) between frames 200 and 
210 is 2.19 in. and this is known to be actually 5.33 ft. 
on the model, the average speed between frames is 


Vir = 5.33 ft./2.19 in. X 3.35 in./0.417 sec. = 19.5 


This speed is plotted for the average frame 205. 

While the entire theory of the method used is covered 
here, it should be pointed out that, if the model path 
could be maintained at 90° to the camera axis at all 
times, there would be no necessity for the foreshortened 
correction and trim could be read directly. Actually, 
it has been found that the angle can vary as much as 
10° without introducing more than 0.5° error in trim. 
With increased experience in handling the model, it has 
been found possible to maintain the course +10° to 
the camera axis, which has greatly simplified the anal- 
ysis. The timesaving through maintaining course 
to these tolerances justifies rerunning a run if the toler- 
ances are exceeded. This procedure has now become 
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XP4Y-1 RADIO MODEL 
ROUGH WATER TESTS 


RUN NO. 3C 3-16-45 
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Fic, 15, Typical run plot, 


standard practice and the analysis consists only of read- 
ing trim, longitudinal and speed reference. The sample 
computation and plot of Fig. 15, however, includes the 
entire procedure for the record. 

An additional interesting observation of the plot of 
Fig. 15 is that through the use of the angle of y, direc- 
tional stability can be readily determined and plotted 
directly as yaw path versus time or speed. Also, it is 
obvious that accelerations can be determined rapidly, 
without the use of accelerometers, by taking the second 
derivative of the trim versus time plot for rotational 
acceleration, and displacement versus time for transla- 
tional accelerations.‘ Vertical displacement relative 
to the horizon can be determined easily, if desired, in 
the analyzer. , 

The system of instrumentation and analysis herein 
described is convenient for free body analysis, is 
technically sound, and provides a method of analysis 
approaching the wind tunnel in accuracy. The method 
and procedures are simple and direct and do not re- 
quire highly skilled technicians. The versatility of 
the method is obvious. 


CORRELATION STUDIES 


The primary purpose of any new research technique 
or procedure is to establish correlation with some known 
or accepted procedure of a similar nature. This was 


recognized when the original program of radio remote 
control of dynamically similar models was organized 
and the Consolidated Vultee Model 31 airplane was 
chosen as the experimental model. That this was a 
wise choice is borne out by the excellent correlation 
that was obtained, not only with a similar procedure, 
such as the N.A.C.A. Towing Basin, but directly with 
the full-scale airplane. It was fortunate that in 1943 
the N.A.C.A. Events Recorder was installed in the 
Model 31 airplane and a limited set of full-scale stability 
limits were run. The Events Recorder recorded si- 
multaneous readings of flap and elevator position, pro- 
peller torque and propeller r.p.m., air speed, water 
speed, and trim. These data were recorded continu- 
ously on sensitized paper and the analysis and pre- 
liminary plots were made under the supervision of Mr. 
Ebert, of the N.A.C.A. Hydrodynamics Staff (refer- 
ence 5). These preliminary, unpublished plots were 
used directly in the correlation reported herein. 


Reference 1 is an unpublished report of the N.A.C.A. 
Towing Basin on dynamic stability tests of a one- 
eighth-scale powered model of the Model 31 airplane 
similar to the -one-eighth-scale radio model reported 
herein. The trim tracks of reference 1 have been used 
in the subject correlation plots. While the radio 
model incorporated a single tail, wind-tunnel tests con- 
ducted at Massachusetts Institute of Technology indi- 
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cated that the moment characteristics were comparable 
to the twin-tail arrangement, there being a slight in- 
crease in aerodynamic stability and elevator effective- 
ness. This slight difference is thought to be well 
within the experimental accuracy of any of the three 
methods of determining hydrodynamic stability. It 
was not possible to find directly comparable runs in all 
cases; however, a sufficient number of such runs are 
presented to justify the conclusion that the three 
methods are comparable. As the full-scale data are 
for a gross weight of 51,000 lbs., this weight has been 
used throughout for the comparisons. 

Figs. 16 through 20 show plots of comparative trim 
tracks for various c.g. positions bracketing both the 
forward and aft c.g. limits for stability. It can be 
readily noted that the character and amplitude of por- 
poising of the radio model and full scale are similar. 
While the character of the porpoising cycle is not de- 
termined in the towing basin, the amplitude and speed 
range are seen to be in good agreement. 

It should be noted that the trim tracks below the 
hump show a consistent variation. It is not immedi- 
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ately apparent why the towing basin shows consist- 
ently higher trims at a given speed unless this is be- 
cause of the partial restraint of the pivot, which forces 
the model to rotate about the center of gravity, and 
the fact that the thrust between the propellers and the 
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Fic. 21. 


towing carriage is not exactly matched. This is not 
considered serious and does not explain why the radio 
model runs at a consistently lower trim than full scale. 
It is noted that the actual hump trim of the radio model 
is in good agreement with full-scale results. 

Regardless of the higher hump trim shown by the 
towing basin model, it is of interest that the speed and 
trim at which porpoising starts is in close agreement in 
all cases. This is considered extremely important and 
of primary interest. Another interesting character- 
istic, where the restraint of the towing pivot in the 














Correlation of spray patterns. 


towing basin may be masking the true results, is the 
fact that the towing basin tests usually exhibit strong 
damping tendencies after a porpoising cycle is estab- 
lished. This results in the towing-basin model ex- 
hibiting stability at the higher speeds; whereas the 
free bodies, both full-scale and the radio model, show a 
tendency to continue oscillating to getaway. This 
may possibly be the result of the reinforcement of the 
porpoising cycle by the more turbulent water en- 
countered both full-scale and with the radio model. 
In marginal cases of instability, tnis turbulent action 
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is possibly sufficient to prevent complete recovery, as 
is the case in the smooth water of the towing basin. 
This condition has been observed in the towing basin 
if runs are repeated before the waves completely sub- 
side. 

Spray correlation between the radio model and full 
scale has been excellent. Fig. 21 presents some com- 
parative pictures that markedly bring out this fact. 
The radio-model method of testing lends itself admir- 
ably to spray studies as the model can be maneuvered 
into any position for direct study. One point of par- 
ticular importance is the fact that with open water tur- 
bulence present, the character of the spray is not 
masked by smooth water blisters resulting in the form 
of spray directly associated with full scale. An ad- 
vantage, too, is the fact that spray conditions in various 
degrees of water surface turbulence can be readily 
studied, including actual rough water and quartering 
wave systems. Taking all factors into .consideration, 
it is felt that the correlation studies presented com- 
pletely justify the principle of free-body testing. 

While the studies discussed herein have applied to 
a fixed ground control station and apply principally to 
the dynamic stability of seaplanes during take-off and 
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landing, their application to other problems of dyna- 
mics is apparent. From a technical standpoint it is 
entirely feasible to operate the transmitter while air- 
borne, thus extending the altitude and range of opera- 
tion. These and additional applications are currently 
under development. With the full development of 
this procedure it will be possible to determine the dy- 
namics of a new design with the same ease and ac- 
curacy as the statics are now determined in a wind 
tunnel. 
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A Method of Testing Thin Sheet Material 


in Compression 


BAXTER C. MADDEN, JR.+ 
Air Materiel Command 


SUMMARY 


A method of determining the compressive mechanical proper- 
ties of thin sheet material is presented. The test jig, specimen, 
and details of the testing procedure are described, and the detec- 
tion of buckling and the measurement of friction between the 
specimen and the jig are discussed. Asa rule, this friction is low 
enough to be ignored but may be readily determined for each 
specimen. Results of tests by the ‘“‘pack’’ method, the “‘cylia- 
der’’ method, and the method described, on 0.035-in. thick steel, 
are given. These indicated good agreement among the three 
methods. Tests on 0.040-in. thick R301-T aluminum alloy 
show consistent results from specimen to specimen. Results 
of tests of 0.020-in. thick clad 24S-T and 0.016-in. thick clad 
75S-T aluminum alloys are also given. Test results are pre- 
sented as stress-strain diagrams. Ease of specimen preparation, 
simplicity of the jig, and the assurance of low friction between 
the specimen and its supporting plates are advantages of the test 
method described. 


INTRODUCTION 


6 dus DETERMINATION of the compressive mechanical 
properties of sheet material, ranging in thickness 
from 0.020 to 0.250 in., is frequently required in con- 
nection with aircraft materials testing. The “pack” 
method! has been used for a number of years and has 
apparently given satisfactory results; however, the 
time required for the preparation of specimens and 
for assembly in the jig is a serious disadvantage. The 
method herein described permits the relatively rapid 
determination of compressive properties of sheet mate- 
rial and is now employed in the Materials Laboratory 
at Wright Field. 


DESCRIPTION OF JIG AND SPECIMEN 


Most of the components of the jig are shown, par- 
tially assembled, in Fig. 1. The jig consists essen- 
tially of two identical rectangular steel plate members, 
two 5/,.-in. diameter dowel pins, four '/2-in. diameter 
bolts with nuts and washers, and two spacers, appear- 
ing on each vertical margin of the left-hand plate 
member. The flat rectangular specimen is seen in its 
proper position between the two spacers. Also shown 
are two Tuckerman optical strain gages employing 
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0.4-in. lozenges. Fig. 2 shows the surfaces of the plate 
members opposite to those shown in Fig. 1. Illustrated 
in Fig. 3 is the jig completely assembled with the speci- 
men and Tuckerman strain gages in place at the begin- 
ning of a test; only one autocollimator is shown in posi- 
tion, the one normally nearer the reader having been 
removed to present greater detail. An arrangement 
using the jig and Huggenberger strain gages is shown 
in Fig. 4. : 

The rectangular plate members were of mild steel 
and were initially */, in. thick, 4 in. wide, and 5'/, in. 
high. The two dowel-pin holes were drilled and reamed 
and serve to maintain the proper positioning of the 
two plates with respect to each other. All holes were 
laid out symmetrically about both vertical and hori- 
zontal centerlines of the plates, the '/2-in. diameter 





Specimen and components of jig partially assembled 
shown with Tuckerman strain gages. 





Fic. 2. Plate member surfaces opposite those shown in Fig. 1. 
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Fic. 3. Jig and specimen assembled prior to test, using Tucker- 
man strain gages. 


hole centers being apart 2*/s in. vertically and 2'/, 
in. horizontally and the °/i,-in. diameter dowel-pin 
hole centers beirg on the horizontal center line and 
2'/, in. apart. A centrally located opening 15/1 in. 
wide and 1*!/z in. high was machined in each plate. 
Surrounding this opening a recess was machined in the 
surfaces shown in Fig. 1 to a depth of 0.25 in. to form a 
shoulder 1/s in. wide against which the aperture insert 
would set. 

The aperture inserts were of mild steel '/, by 1°/is by 
1*/s. in. with two apertures in each. The apertures 
were 0.12 by 0.80 in., symmetrically located in the 
insert with their horizontal centerlines 1 in. apart. 
As shown in Fig. 2, there is a bevel surrounding the 
apertures, the upper aperture having a 30° bevel and 
the lower aperture having a 45° bevel. The insert 
was fitted against the shoulder and electrically welded 
in place. Following this, the surface to be adjacent to 
the specimen was ground to a smooth plane on each 
plate. The plates were then dowel-pinned and bolted 
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Fic. 4. Jig and specimen assembled prior to test, using Huggen- 
berger strain gages. 


together aud the top and bottom of the assembly 
ground parallel. A portion between the apertures of 
each insert was milled to a depth of approximately 
0.030 in., as shown in Fig. 2, to provide clearance for 
the body of the Tuckerman strain gage. 

The four adjustable clips, shown in Fig. 2, each fast- 
ened by a screw and washer, are merely to anchor an 
elastic band for holding the Tuckerman gages in place, 
when these gages are used, and their dimensions are 
such as to allow sufficient clearance around the gages. 

The specimen is a rectangle, with machined edges, 
1.80 in. wide and a length equal to the height of the 
jig plus 0.0€0 to 0.080 in. If the specimen blank is to 
be sheared, an allowance at least equal to the sheet 
thickness is made on each edge for machining to the 
finished dimensions in order to eliminate any cold-work. 
effect due to shearing. Any burrs after machining are 
carefully removed without chamfering the edges. 


TESTING PROCEDURE 


The thickness of the specimen is determined at a 
number of locations, the average being used in deter- 
mining the cross-sectional area. For specimens 0.020 
to 0.040 in. thick, the shims are adjusted to a thickness 
of 0.0005 to 0.0015 iv. greater than the maximum 
measured thickness of the specimen; for specimens 
0.040 to 0.064 in. thick, the total clearance should be 
0.0015 to 0.0030 in.; for specimens 0.064 to 0.125 in. 
thick, the total clearance should be 0.003 to 0.004 in; 
for specimens 0.125 to 0.250 in. thick, the clearance 
should be 0.004 to 0.005 in. The shim thickness is 
conveniently obtained by using the same material 
as that being tested and by placing one or more layers 
of paper over each metal shim. The paper laminations 
are quickly made by placing them in position between 
the two plates and punching out the holes with sec- 
tions of '/2- and 5/,.-in. diameter rod, the end of each 
being machined slightly concave. 

Both sides of the specimen are thoroughly coated 
with lead soap base lubricating grease (U.S. Army 
Specification 2-118). This appears to be most easily 
applied by smearing on with the fingers. It is im- 
portant that the grease and the surfaces of the specimen 
and jig be kept free of any grit or other foreign particles. 
Various oils, waxes, and other greases were tried as 
lubricants and found to be unsatisfactory because of 
insufficient film strength or excessive friction. The 
jig is assembled with the specimen in position as shown 
in Fig. 1. It is convenient to lay one plate with the 
bolts and dowel pins in place on a horizontal surface. 
The shims and specimen are then positioned. Finally, 
the other plate is fitted to the bolts and dowel pins, 
the washers are installed, and the nuts are evenly tight- 
ened. 

The specimen is moved back and forth two or three 
times so that opposite ends of the specimen alternately 
become flush with the ends of the jig. The jig is then 
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placed upright on supports with a vacant space directly 
beneath the specimen and the bottom of the specimen 
flush with the bottom of the jig. Small weights are 
carefully balanced on the top of the specimen. A load 
of 5 lbs. should overcome the friction due to the grease 
and cause the specimen to sink slowly until its top is 
even with that of the jig. If the load required is greater 
than 5 lbs., this indicates that a burr remains on the 
specimen, that a foreign particle is between the jig and 
the specimen, or that insufficient clearance exists be- 
tween the surfaces of the jig and the specimen. Such 
condition should be corrected before proceeding with the 
test. Having determined that excessive friction is not 
present, the specimen is set with its lower edge flush 
with the bottom of the jig. 

The strain gages are installed with the specimen 
contacts against the upper edge of each aperture or 
nearly so. This ensures that, as downward deforma- 
tion of the specimen occurs, the specimen contacts 
will not be hindered in their travel by the lower edges of 
the apertures. When Tuckerman gages are used, the 
movable contact is placed in the upper aperture where 
the bevel is 30°; this allows the clearance necessary 
for the movement of the lozenge. 

The assembled jig with specimen is placed upon a 
ground steel plate at the center of the lower platen of 
the testing machine, as shown in Fig. 3. A piece of 
paper between the bottom of the jig and the plate 
serves to minimize any uneven loading that may re- 
sult from the end of the specimen not being exactly 
true or from minute machining marks on the end of the 
specimen. A piece of metal, approximately 0.040 by 1 
by 2 in. and somewhat softer than that being tested, is 
approximately centered on the top of the specimen, as 
indicated in Fig. 3, just beneath the spherical seat 
arrangement, which is centered with respect to the 
top of the specimen. For relatively thin specimens, 
an aid in centering the spherical seat arrangement 
is to balance it on the top edge of the specimen. Initial 
readingsof the strain gages are taken. A load of approxi- 
mately 2 per cent of the expected yield strength of the 
specimen is applied and strain-gage readings are again 
taken. Ifthe strain increments on each side of the speci- 
men are approximately equal, the test is continued; if 
not, the load is removed and the spherical seat is moved 
slightly toward the side of the specimen having the 
lesser strain increment. The test for equal strains on 
each side of the specimen is repeated and further adjust- 
ment made if necessary. The load-strain data are then 
obtained by reading simultaneously both strain gages 
and the load indicator as the load is slowly increased to 
somewhat beyond the yield strength of the material. 


DETECTION OF BUCKLING OF THE SPECIMEN 


The determination of the compressive yield strength 
on sheet material is valid only if the specimen surface 
remains free of buckles. Should buckling occur, its 


commencement may be readily detected by comparing 
the increments in strain-gage readings. While the 
specimen surface is plane and the specimen is under an 
axial load, a consideration of pairs of simultaneous 
strain-gage readings shows that the increments indi- 
cated by one gage are nearly a constant proportion of 
the corresponding increments indicated by the other 
gage, except for an occasional discrepancy probably 
resulting from instability of the testing machine; differ- 
ences in gage factors account for the corresponding in- 
crements generally not being approximately equal. As 
a buckle forms, this proportionality is destroyed and 
remains so for subsequent readings. The cause of this 
may be explained by reference to Fig. 5, where the com- 
mencement of a buckle is shown exaggerated. 

Whenever a buckle has developed, it has been ob- 
served to occur at one of the apertures because of lack 
of support. Therefore, the early stage of a buckle is 
shown at an aperture, where one contact of each 
strain gage bears on the specimen. The distance be- 
tween the contacts of the strain gage on the convex 
side of the buckle is indicated as L1, while that for the 
gage on the concave side is indicated as L2. Prior 
to the beginning of a buckle, L1 and L2 are substanti- 
ally equal for any given load. The effect of the buckle, 
taken alone, is to produce a compressive strain on the 
concave surface and a tensile strain on the convex sur- 
face. 

The net result of this, when considered along 
with the normal compressive strain, is that more axial 





STRAIN GAGE CONTACT 








/ SPECIMEN 

















“ir 
STRAIN GAGE CONTACT 


Fic. 5. Commencement of a buckle in a specimen. 
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deformation exists on the concave surface than on the 
convex surface and that L2 becomes substantially less 
than L1, with the effect continuing to become more ap- 
parent as the size of the buckle increases. The original 
data in Table 1 from the compression test of a specimen 
0.009 in. thick (16 per cent manganese, 12 per cent 
chromium steel), where buckling occurred, illustrate the 
behavior of the strain gages. Gage No. 268 was in 
contact with the concave side of the buckle. At loads 
above 1,550 Ibs., buckling was evidenced by the per- 
sistent increase in the increment ratio. 








TABLE 1 
—Strain Gage No..268—. —Strain Gage No. 378— 
Incre- Incre- 
ment ment Increment 
’ over over Ratio 
Load, Gage previous Gage previous (No. 268+ 
Lbs. reading reading reading reading No. 378) 
100 0 ta 0 — 
150 22 22 26 26 0.85 
200 52 30 56 30 1.00 
250 82 30 84 28 1.07 
300 112 30 112 28 1.07 
850 140 28 138 26 1.08 
400 174 34 172 34 1.00 
450 204 30 202 30 1.00 
500 236 32 234 32 1.00 
550 268 32 266 32 1.00 
600 304 3 300 34 1.06 
650 334 30 330 30 1.00 
700 368 34 362 32 1.06 
750 402 34 396 34 1.00 
800 436 34 426 30 1.13 
850 472 36 462 36 1.00 
900 506 34 496 * 34 1.00 
950 544 38 532 36 1.06 
1,000 582 38 562 30 1.27 
1,050 620 38 598 36 1.06 
1,100 660 40 640 42 0.95 
1,150 700 40 680 40 1.00 
1,200 740 40 720 40 1.00 
1,250 780 40 760 40 1.00 
1,300 820 40 800 40 1.00 
1,350 870 50 840 40 1,25 
1,400 930 60 900 60 1.00 
1,450 990 60 960 60 1.00 
1,500 1,040 50 1,020 60 0.83 
1,550 1,100 60 1,080 60 1.00 
1,600 1,170 70 1,140 60 1.17 
1,650 1,260 90 1,220 80 1.13 
1,700 1,370 110 1,300 80 1.38 
1,750 1,530 160 1,460 100 1.60 
1,800 1,720 190 1,520 120 1.58 
1,850 1,960 240 1,680 160 1. 50 


A definite check as to whether buckling occurred is 
had when the jig is disassembled after a test. With no 
buckling, the surface of the specimen appears as shown 
in Fig. 6, having no breaks in the grease coating. If 
even slight buckling has occurred, the specimen surface 
appears as shown in Fig. 7, with areas where the grease 
film has failed and has allowed metal to metal contact. 
A severe buckle appears as in Fig. 8, where two views 
of a buckled specimen are shown. 


DETERMINATION OF FRICTION BETWEEN SPECIMEN AND 


Jic 


In any compression test where the specimen must be 
supported to increase its stability against buckling, 





Fic. 6. Specimen surface after test. No buckling occurred 





Fic. 7. Specimen surface after test. Slight buckling occurred. 
friction due to such support becomes a factor that 
should be evaluated. 

As the specimen is compressed, there occurs an in- 
crease in thickness which may be calculated from Pois- 


son’s ratio. An average of twelve determinations? of 
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Fic. 8. Two views of a specimen after test. Severe buckling 
occurred. 
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Fic. 9. Cross section showing arrangement for measuring 
maximum friction load during compression test. 


Poisson’s ratio at various strains gave a value of 0.352 at 
the 0.2 per cent offset for various aluminum-alloy 
sheet materials. For an aluminum alloy with a 0.2 
per cent offset compressive yield strength of 70,000 and 
10,500,000 Ibs. per sq.in. modulus of elasticity, the 
total strain at the yield strength is (70,000/10,500,000) 
+ 0.002 = 0.0086 in. per in.; the increase in thickness 
= (0.352) (0.0086) = 0.00303 in. per in. This re- 
sults in approximate increases in thickness of 0.00006 in. 
for 0.020-in. thick material, 0.0002 in. for 0.064-in. 
thick material, and 0.0008 in. for 0.250-in. thick mate- 
rial. Approximately the same increases in thickness 
should occur with magnesium alloys and steels when 
stressed to their respective yield strengths. 

The friction, which is determined not to exceed 5 
Ibs. after tightening the bolts and nuts on the jig, 


would merely ensure that the technique up to that 
point had been correct. The friction during the test 
may be expected to increase somewhat because of the 
increase in thickness noted above. 

To determine the maximum friction during the test, 
the arrangement shown in Fig. 9 was employed. The 
weight of the jig plus the friction load was transmitted 
through two !/,-in. thick 14S-T aluminum-alloy plates, 
each in contact with a 1-mm. diameter steel ball rest- 
ing in a drilled center mark in a '/s-in. thick steel strip, 
The jig was carefully set so that the ball centers were 
in the plane bisecting the thickness of the specimen, 
The test was then run to a load corresponding to a 
stress somewhat above the yield strength of the mate- 
rial. After completing the compression test, the diame- 
ter of the impression in each aluminum plate, due to 
the steel ball, was measured with the measuring micro- 
scope on a Vickers-Armstrongs hardness tester. Using 
this machine, a series of hardness tests, with a 1-mm, 
diameter steel ball and various loads, was made on the 
aluminum-alloy plates, resulting in the plotting of a 
curve relating diameter of ball impression to load on the 
ball. The load taken by each ball during the compres- 
sion test was determined from this curve. The sum 
of the two ball loads, minus the weight of the jig and 
aluminum plates, represented the maximum fric- 
tion load. The following data were obtained in this 
test: 


Specimen: R301-T aluminum alloy, with long axis perpendicu- 
lar to rolling direction, 0.0619 in. thick, 1.8004 in. wide, 63,550 
Ibs. per sq.in. compressive yield strength at 0.2 per cent 
offset. 

Maximum load on specimen during test was 7,200 Ibs., corres- 
ponding to 64,600 lbs. per sq.in. 

Total clearance between surfaces of specimen and jig—0.0025 
in. 

Friction load prior to test—3.0 Ibs. 

Weight of jig and aluminum plates—8.9 Ibs. 

Maximum total load on two 1-mm. diameter balls during com- 
pression test—23.4 lbs. 

Maximum friction load during compression test—14.5 Ibs 

Room temperature during test—74°F. 


It appears reasonable to assume that the friction 
load is distributed uniformly over the two surfaces of the 
specimen. The lower half of the specimen may be 
considered as contributing, through friction, one-half 
the friction load; this is effective iv producing strain 
at the strain gages. The upper half of the specimen 
may be considered to contribute through friction one- 
half the friction load; however, this is not effective in 
producing strain at the strain gages. Therefore, the 
maximum load on the specimen, as indicated by the 
testing machine, may be considered as being too high 
by an amount equal to one-half the maximum friction 
load occurring during the test. In the test cited above, 
this error would be 7.3 lbs., which would usually be 
negligible. However, should it be desirable to evaluate 
the friction load on each test, this may be done with- 
out difficulty. 





lel 
loa: 
Tw 
spe 
stra 
twc 
The 
yiel 
tair 


pres 
allo 
wit] 
and 
the 

fron 
chin 
shor 


alloy 
sent 


thic] 


the | 


critic 
of 0. 
buck 
perp 
buck 
to tl 


sq.in 


used 
mini: 
press 
ness) 


that 
> test 
f the 


- test, 
The 
Litted 
lates, 
rest- 
strip. 
were 
men, 
toa 
nate- 
ame- 
ue to 
licro- 
Jsing 
-mm, 
n the 
of a 
n the 
pres- 
sum 
and 
fric- 
this 


idicu- 
3,550 
cent 


prres- 


.0025 


com- 


‘tion 


f the 








TESTING THIN SHEET MATERIAL IN COMPRESSION 351 


COMPARISON WITH ‘“PACK’’ METHOD AND ‘‘CYLINDER”’ 
METHOD 


A comparison of the method described with two 
other methods of determining compressive mechanical 
properties of sheet meterial was made. A strip of 
hard-rolled 16 per cent manganese steel, 0.035 in. thick, 
was made into a number of specimens for test by the 
“pack’’ method,’ the “cylinder” method,* and the 
method herein described. One-half of the specimens 
were prepared so that the load would be applied paral- 
lel to the direction of rolling; the remainder, so that the 
load would be applied perpendicular to this direction. 
Two Tuckerman optical strain gages were used on each 
specimen during the test. In Fig. 10, typical stress- 
strain diagrams are presented for tests in each of the 
two directions, as obtained by the three methods. 
These diagrams indicate that the method described 
yields results which compare very well with those ob- 
tained by either of the other methods. 


RESULTS OF TESTS ON VARIOUS ALUMINUM ALLOYS 


To check upon the accuracy of duplication of com- 
pressive stress-strain diagrams with a given aluminum 
alloy and by the test method described, two specimens 
with the long axis parallel to the direction of rolling 
and two specimens with the long axis perpendicular to 
the direction of rolling, cut adjacent to one another 
from a sheet of 0.040-in. thick R301-T alloy, were ma- 
chined and tested. The results were plotted and are 
shown in Fig. 11. It is evident that for each pair of 
curves the second is nearly a duplicate of the first. 


Compressive stress-strain diagrams, determined by 
the method described, for 0.020-in. thick clad 24S-T 
alloy and for 0.016-in. thick clad 75S-T alloy are pre- 
sented in Fig. 12. It is to be noted that the 0.016-in. 
thick clad 75S-T alloy, loaded perpendicular to the 
direction of rolling, failed because of buckling before 
the 0.2 per cent offset was reached. 


LIMITATIONS ON THICKNESS OF SPECIMENS 


The tests of the clad 75S-T alloy denote that the 
critical buckling stress for this material in a thickness 
of 0.016 in. is approximately 70,000 Ibs. per sq.in., since 
buckling occurred at this stress in the specimen loaded 
perpendicular to the rolling direction and since no 
buckling occurred in the specimen loaded parallel 
to the rolling direction at a stress of 69,200 Ibs. per 
sq.in. 

For each of the three general types of alloys widely 
used in thin gages in aircraft, au indication of the 
minimum thickness that may be loaded to the com- 
pressive yield strength without buckling (critical thick- 
ness) may be had by use of the following relation- 
ship :* 

F, = KE(t/b)? (1) 
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where 
F., = critical stress at which buckling occurs 
at one of the apertures 
E = modulus of elasticity 
t = thickness of specimen 
b = length of aperture perpendicular to direction of 


load 
K = a factor depending on the type of loading, 
dimensions of aperture, edge fixity condi- 
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tions at the aperture, and Poisson’s ratio. 
Since these items are substantially constant 
for the three types of alloys, K may be ap- 
proximately evaluated from the test on 
0.016-in. thick clad 75S-T aluminum 
alloy. 


Then 


F. 70,000 


ae ae - — = 16.7 (2) 
E(t/b)? 10,500,000 (0.016/0.80)? 





The critical thickness, from Eq. (1), if F,, is taken as 
the yield strength of the material, becomes 


t = bV F,,/KE = 0.196V F,/E (3) 


The highest approximate compressive yield 
strength (at 0.2 per cent offset) likely to be en- 
countered at present in each type of alloy is shown in 
Table 2. 





TABLE 2 





Compressive Yield Modulus of Elas- 
Strength, Lbs. per ticity, Lbs. per 


General Type of Alloy Sq.In. $q.In. 

Corrosion resisting steel 200,000 28,000,000 
Aluminum 80,000 10,500,000 
Magnesium 40,000 6,500,000 


By substituting the above values for /,,and £ in Eq. 
(3), the following critical thickness of specimen for 
each type of alloy is determined: corrosion resisting 
steels—0.017 in.; aluminum alloys—0.017 in.; and 
magnesium alloys—0.015 in. . 

For a material whose compressive yield strength is 
less than that assumed -for the general type of alloy, 
thinner critical gages would be indicated. It is to be 
understood that Eq. (1) is predicated upon stresses not 
exceeding the proportional limit; therefore, the above 
calculations are rather approximate. However, a lower 
limit of specimen thickness of 0.020 in. appears to be 
conservative. 


ADVANTAGES OF METHOD 


Compared with the “‘pack’’! and ‘‘cylinder’’* methods, 
the ease of specimen preparation prior to the test ap- 
pears to be one of the important advantages of any of 
the methods of compression testing of thin sheet mate- 
rial wherein a single thickness of the flat sheet is used as 
a specimen. 

A further advantage of the method described over 
other single thickness methods, wherein supporting 
plates bear against two sides of the specimen and the 
bearing force is applied by a screw, set merely accord- 
ing to the ‘‘feel’’ of the operator, is that a definite 
amount of clearance is allowed for the lubricant film, 
resulting in low friction between the specimen and sup- 
porting plates during the test. It has been shown! 
that considerable error in both yield strength and modu- 
lus of elasticity may occur in compression tests in 
which the manually adjusted initial supporting force is 
too large; both of these properties appear to increase 
as the initial supporting force is increased. 

The, jig used with the test method described appears 
to be considerably simpler to construct and maintain 
than the fixture used where a flat, single thickness 
specimen is supported by two series of rollers.* 
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Helicopter Performance 


ROBERT L. LICHTEN* 
Kellett Aircraft Corporation 


SUMMARY 


A straightforward method for calculating helicopter perform- 
ance is developed which retains the accuracy of the blade-element 
approach yet is not prohibitively lengthy. 

Expressions are presented for the flapping, thrust, longitudinal- 
force, and power coefficients of a hinged-blade lifting rotor, which 
include the effect of linearly tapered and twisted blades. These 
expressions have been simplified for use by presenting them as 
simple functions of the inflow ratio and the blade pitch angle, 
having numerical coefficients. Tables of the necessary numerical 
coefficients which are satisfactory for most applications are pre- 
sented in the Appendix. 

Performance is expressed in terms of the climb ratio, an ex- 
pression for which is developed based on the force equilibrium 
existing in steady flight. The climb ratio may be presented for a 
given helicopter in a series of charts covering the speed range. 
These charts will express the complete performance of the heli- 
copter. A cross-plot from these charts yields a single chart 
which expresses the complete level flight performance. 

Included in the paper are examples of such charts, a compari- 
son between the calculated and flight-test performance of a heli- 
copter, and curves showing the effects on performance of various 
changes in an example helicopter. 


INTRODUCTION 


: he AERODYNAMIC ANALYSIS of the lifting rotor in 
forward flight by integration of the forces acting on 
a blade element has been developed by many authors.'~* 
As successive analyses became more and more exact, 
they tended to become more and more complex, so that 
performance calculations for rotary wing aircraft using 
these analyses became quite lengthy and difficult. As 
a result, there has been a tendency to calculate per- 
formance by approximate methods.?—"* Such methods 
have generally. utilized greatly simplified forms of the 
blade-element expressions for rotor characteristics, 
adapted and combined to yield the total rate of energy 
consumption for the helicopter in flight. 

When applied by experienced persons, satisfactory 
performance evaluations may be made by approxi- 
mate methods. However, performance calculated by 
such methods is likely to be optimistic in the cruising 
and high-speed range. In addition, some approximate 
methods do not permit determination of the effects of 
changes in blade pitch, rotor speed, rotor solidity, and 
blade plan form, while others are not accurate for cal- 
culation of climb performance. 

A notable advance in simplifying calculations using 
the more exact blade-element expressions was made by 


Presented at the Philadelphia Section of the I.A.S., Phila- 
delphia, December 18, 1945. Received April 5, 1946. 
* Aerodynamicist. 


Bailey. He introduced the use of numerical coeffi- 
cients, which are valid for most rotors of conventional 
design, in the calculation of rotor aerodynamic char- 
acteristics, thus greatly simplifying numerical applica- 
tion of the expressions. Bailey’s analysis is limited to 
the case of untapered rotor blades. More recent pub- 
lications of the N.A.C.A. have suggested a method for 
the use of his analysis in calculating helicopter per- 
formance. The method proposed is somewhat un- 
wieldy, the speed range being only partially covered, 
and a method of successive approximations being re- 
quired to select desired values from a set of calculated 
charts. 

In the present paper, expressions for rotor character- 
istics are presented which include the effect of blade 
taper as well as twist. Numerical coefficients of the 
type introduced by Bailey are tabulated in Tables 1-4 
for use in evaluation of the rotor expressions. The drag 
of the rotor is treated by means of the longitudinal force, 
H. 

‘Thrust and longitudinal-force coefficients together 
with the other parameters of the helicopter are com- 
bined into a simple expression for the climb ratio x, 
which is a fundamental variable of performance. The 
climb ratio is defined as the ratio of rate of climb to tip 
speed. Values of the climb ratio and rotor power coef- 
ficient (identical numerically to torque coefficient) may 
be computed for the range of possible operating condi- 
tions for a particular helicopter. These values may be 
presented in a series of charts, examples of which are 
given later in the paper. From such charts, desired 
performance items for the helicopter may be deter- 
mined rapidly and easily. 

The general method used in the derivation of the 
rotor expressions is the same as that presented by 
Wheatley ;* therefore, the derivations are not given in 
detail in this paper. One difference is that second har- 
monic flapping was assumed to be negligible for rotors 
having blades of present-practice mass characteristics, 
insofar as the effect of second harmonic flapping on the 
thrust, torque, and longitudinal force of the rotor is 
concerned. 

The analysis is valid for helicopters having flapping- 
or teetering-type rotors which transmit little or no 
moment to the rotor hub. It may also be adapted for 
use with cantilever-type rotors, provided that such 
rotors are operated with hub moments substantially 
reduced to zero in all flight conditions by means of cyclic 
feathering. 
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DEFINITIONS 


The virtual axis of rotation is an axis about which the 
rotor hub would have to rotate in order to maintain 
unchanged the path of the rotor blades in space, for 
any given steady flight condition, if all cyclic feathering 
control were eliminated from the rotor. 

The plane of the rotor disc is a plane perpendicular to 
the virtual axis of rotation, passing through the rotor 
hub. 


NOTATION 


= rotor disc area = wR? for single rotor, sq.ft. 


4 
— 
| 


a = slope of the blade section lift coefficient curve, per rad. 
a, = coefficient in Fourier series for blade flapping motion 
B- = tip-loss factor 


b = number of blades per retor 

6, = coefficient in Fourier series for blade flapping motion 
£ = blade chord at any radius = ¢o(1 + #x), ft. 

projected blade chord at center of rotor, ft. 


Lo = 
Cp, = helicopter parasite drag coefficient = f/A 

C, = blade section lift coefficient = aa, 

Cx = rotor longitudinal-force coefficient = H/[!/2p(wR)?A | 


Cp = rotor shaft power coefficient = P/[!/29(wR)#A | 

Cr = rotor thrust coefficient = 7/[!/:p(wR)?A } 

D; = parasite drag of all parts of helicopter except main 
rotor blades, Ibs. 

f = helicopter parasite drag area = Dy/ (2/29 V?), sq.ft. 

F, = numerical coefficient in simplified expressions for flap- 
ping coefficients 

H = rotor longitudinal force, vector lying in the plane of 
the rotor disc directed along the azimuth y = 0, Ibs. 

H, = numerical coefficient in simplified expression for Cy 

I, = moment of inertia of one blade about its flapping axis, 
slug-ft.? 

K = parameter defined by Eq. (3) 

kg = rate of change of blade pitch with flapping angle = 


d6/dg 

Mw = weight moment of one blade about its flapping axis, 
ft.lbs. 

P = rotor shaft power, ft.lbs. per sec. 

P, = numerical coefficient in simplified expression for Cp 

r = radius toa blade element, ft. 

R_ = radius to blade tip, ft. 

t = blade taper factor = (cyp/co) — 1 

T = rotor thrust force, vector coincident with virtual axis 
of rotation, directed upward, lbs. 

Tx = numerical coefficient in simplified expression for C7 

V = true air speed along flight path, ft. per sec. 

V’ = resultant velocity at rotor disc, ft. per sec. 

v = induced velocity at rotor disc directed opposite to 
thrust vector, ft. per sec. 

v; = additional induced velocity varying linearly from front 
‘to rear of rotor disc, ft. per sec. 

W = gross weight of helicopter, lbs. 

x = nondimensional radius to a blade element = r/R 

a = rotor angle of attack, the angle between the flight path 


and the plane of the rotor disc, positive when ad- 
vancing edge of disc is raised, rad. — 

= angle of attack of a blade element, rad. 

B = flapping angle between blade span axis and plane of 
rotor disc, positive upward, rad. 

= blade mass constant =( !/spacoR*) /I; 


5 = blade section profile-drag coefficient = 59 + d:a0 + 
bean? 
$5, = radial blade drag coefficient = twice the coefficient of 


skin friction for flow along blade in radial direction 


0 = pitch angle of a blade element = @ + 6,x, rad. 


‘ = angle of inclination of virtual rotor axis to vertical, 
positive forward, rad. 

K = climb ratio = rate of climb/tip speed = 
—[V sin (a +1)]/wR 

r = inflow ratio = (V sin a — v)/wR 

1 = additional inflow ratio defined by Eq. (3) 

a = advance ratio = (V cos a)/wR 

p = air density, slugs per ft.* 

o9 = rotor solidity based on blade root chord = bceo/xR 

T = rotor blade weight moment parameter = My/I,w? 

y = blade azimuth angle measured from downwind position, 
rad. 

w = rotor rotational speed = dy /dt, rad. per sec. 

ANALYSIS 


Derivation of Rotor Expressions 

The expression for the constant induced velocity 
generated by a single isolated rotor as used by Wheatley‘ 
is 

v = T/2xR*pV’ 

This expression represents an assumption of the magni- 
tude of the induced velocity at the rotor disc, based on 
simple momentum theory and an analogy with finite 
airfoil theory. It has proved to be fairly satisfactory, 
However, it should be noted that for helicopter con- 
figurations other than single rotor, this expression and 
those derived from it may have to be suitably modified. 
In order to be consistent with the use of the tip-loss fac- 
tor B in the rotor expressions, the above expression is 
modified for this paper by use of the effective disc 
radius BR; resulting in 


v = T/2nB?R*pV' (1) 

Accordingly, the expression for the rotor angle of at- 
tack becomes 

tan a = A/u + C,/4B%u(d? + p?)'? (2) 


Recent studies by the N.A.C.A.! and others! have 
tended to confirm the hypothesis of Wheatley‘ that an 
additional induced velocity varying from front to 
rear of the disc should be superposed upon the in- 
duced velocity v. The additional velocity is expressed 


as 

v, = Kvox cos y = \ywRx cos p (3) 
where K is a parameter varying with the flight condi- 
tion from zero at hovering to a probable maximum value 
of unity. For a single rotor, 


dh. = KC,/4B*(? + y?)'* (4) 

The blade pitch angle, including a linear twist from 

root to tip and pitch change with flapping, is given by 

6 = 0 + OxK + kgB (5) 

The blade chord, including a linear taper from root to 
tip, is given by 

c = co(1 + tx) (6) 
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where ¢ is the projected blade chord at the rotor center 
and the taper factor ¢ is defined by 


t = (Ctip/co) — 1 (7) 


The variation with blade-section angle of attack of the 
,section profile-drag coefficient is represented by 


6 = 5 + bia, + boa,” (8) 
The section lift coefficient at a blade element is given by 
CG. = aa (9) 


In accordance with the previously stated assumption 
that second-harmonic flapping may be neglected, the 


assumed form for the Fourier series expressing the blade 


flapping motion is 


B = a — a, cos y — hj sin y (10) 


The next step in the analysis is to determine the 
resultant velocity components at the blade element and, 
from these and the blade section characteristics, the 
aerodynamic forces acting on the blade element. 


to thrust about the flapping axis) over the time interval 
required for a revolution of the rotor. The blade mo- 
tion is determined by solving the differential equation 
of angular motion about the flapping axis. This part 
of the analysis has been extensively treated by other 
authors, who are referred to below, and therefore the 
derivations are not given in this paper. The final ex- 
pressions derived for the flapping, thrust, longitudinal- 
force, and power coefficients are presented in the Ap- - 
pendix. These expressions include the effect of the 
reversed velocity region on the retreating blades as 
developed by Wheatley.‘ 


The simplified forms of the rotor expressions are given 
in the following sections. In order to avoid com- 
plexity, it was found necessary to set kz = 0 for these 
expressions. The effect of this restriction is discussed 
later. 


Flapping Coefficients 
The procedure used to evaluate the flapping coeffi- 































































































These elemental forces are then integrated over the cients is identical to that used by Wheatley.‘ The 
blade span and averaged (except for the moment due simplified expressions are 
TABLE | bee eam rent 
Numerical Coefficients for Equations for Flapping Coefficients 
t 0 —0.4 —0.8 0 —0.4 —0.8 0 —0.4 —0.8 
F, te ne Fy F, 
p=0 0.3042 0.2157 0.1272 0.2213 0.1526 0.0839 0.1718 0.1163 0.0608 
0.1 0.3043 0.2158 0.1273 0. 2237 0.1544 0.0851 0.1733 0.1173 0.0614 
0.2 0.3048 0.2162 0.1277 0.2306 0.1594 0.0883 0.1779 0.1206 0.0633 
0.3 0.3063 0.2177 0.1290 0.2422 0.1680 0.0938 0.1855 0.1260 0.0665 
0.4 0.3093 0.2205 0.1316 0.2581 0.1797 0.1012 0.1961 0.1335 0.0709 . 
0.5 0.3141 0.2248 0.1355 0.2781 0.1942 0.1103 0.2098 0.14382 0.0766 
— F, F; Fs 
n=0 0 0 0 0 0 0 0 0 0 
0.1 0.2135 0.2294 0.2725 0.2762 0.2841 0.3049 0.2013 0.2014 0.2017 
0.2 0.4298 0.4614 0.5453 0.5623 0.5795 0.6250 0.4086 0.4091 0.4105 
0.3 0.6540 0.7008 0.8272 0.8681 0.8970 0.9746 0.6303 0.6330 0.6391 
0.4 0.8899 0.9524 1.1230 1.2089 1. 2545 1.3783 0.8746 0.8803 0.8984 
0.5 1.1449 1.2248 1.4448 1.6024 1.6743 1.8723 1.1532 1. 1667 1.2037 
F, F; seen - — Fy 
p=0 0 0 0 0 0 0 0 0 0 
0.1 0.0416 0.0304 0.0192 0.0306 0.0217 0.0128 0.0237 0.0165 0.0092 
0.2 0.0821 0.0598 0.0378 ° 0.0621 0.0441 0.0261 0.0479 0.0334 0.0187 
0.3 0.1210 0.0882 0.0558 0.0957 0.0680 0.0406 0.0733 0.0510 0.0288 
0.4 0.1579 0.1155 0.0734 0.1318 0.0941 0.0564 0.1001 0.0699 0.0395 
0.5 0.1933 0.1418 0.0907 0.1712 0.1225 0.0739 0.1291 0.0903 0.0513 
Fo = Fy ne 
im 0 0 0. 1.0000 1.0000 1.0000 
0.1 —0.1366 —0.1407 —0.1505 0.9946 0.9941 0.9941 
0.2 —0.2695 —0.2767 —0.2958 0.9792 0.9776 0.9733 
0.3 —0.3950 —0.4050 —0.4326 0.9543 0.9508 0.9427 
0.4 —0.5106 —0.5237 —0.5576 0.9217 0.9165 0.9031 
0.5 —0.6155 —0.6307 —0.6697 0.8827 0.8750 0.8552 
Paley ST ce TABLE 2 ee rear 
Numerical Coefficients for Equation for Thrust Coefficient 
t 0 —0.4 —0.8 0 —0.4 —0.8 0 -0.4 08 
T; Te Ts 
s=0 0.4705 0.3488 0.2271 0.3042 0.2157 0.1272 0.2213 0. 1526 0.0839 
0.1 0.4730 0.3513 0.2296 0.3090 0.2195 0.1300 0. 2237 0.1544 0.0851 
0.2 0.4809 0.3591 0.2373 0.3231 0.2309 0. 1386 0.2310 0.1598 0.0887 
0.3 0.4953 0.3730 0.2610 0.3470 0.2502 0.1536 0.2443 0.1701 0.0960 
0.4 0.5179 0.3950 0.2732 0.3823 0.2795 0.1771 0.2651 0.1868 0.1084 
0.5 0.5521 0.4286 0.3080 ~ 0.4326 0.3224 0.2138 0.2961 0.2124 0.1291 
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a TasLE3 ee ene Sere poy one 7 
Nume rical Coe ficients for Equation for Power Coefficient 
ang 0 —0.4 -~0.8 0 —0.4 —0.8 0 —0.4 0.8 ’ 
a etn cn Py “ —— -P; il 
w= 0 0). 2500 0.1700 0.0900 0.3333 0.2333 0.1333 0.2500 0.1700 0.0900 * 
0.1 0). 2525 0.1718 0.0912 0.3333 0.2333 0.1333 0.2525 0.1718 0.0912 
0.2 0.2599 0.1772 0.0946 0.3333 0. 2333 0. 1333 0.2600 0.1773 0.0947 ° 
0.3 0.2722 0. 1862 0. 1002 0.3333 0. 2333 0.1333 0.2725 0. 1865 0.1005 
0.4 0. 2892 0. 1985 0.1079 0.3333 0. 2333 0. 1333 0.2900 0.1993 0. 1087 
0.5 0.3105 0.2138 0.1172 0.3333 0. 2333 0. 1333 0.3125 0.2158 0.1192 
ak preasisriecsee Siedinen! ” ‘coat arte AAT NE le < Be IM src-taket a _ Ses a 
uw=0 0.2000 0.1333 0.0666 0. 5000 0.3667 ). 2334 0.4705 0.3488 0.2271 m 
0.1 0.2017 0.1345 0.0674 0.5145 0.3776 2411 0.4838 0.3590 0.2347 
0.2 0. 2067 0.1880 0.0693 0.5589 0.4113 0.2645 0.5242 0.3900 0.2571 
0.3 0.2150 0.1488 0.0727 0.6362 0.4697 pabnen 0.5940 ae 0.2948 
0.4 0. 2267 0.1520 0.0773 0.7523 0.5556 0.362% 0.6990 . 5211 0.3488 
0.5 0.2417 0.1625 0. ). 0834 0.9181 0.6766 0. 4408 0.8484 0. 6310 0.4225 
sie BMS oR TOS ren ea estes: “noamasmienican cearneien eel 
u=0 0.6667 0. 4667 0. 2667 0.3042 0.2157 0.1272 0.5000 0.3400 0.1800 » 
0.1 0.6965 0.4892 0.2819 0.3316 0. 2267 0. 1422 0.5217 0.3559 0.1900 
0.2 0.7904 0.5605 0.3312 0.4168 0.3023 0.1886 0. 5894 0.4060 0.2218 
0.3 0.9638 0.6920 0.4217 0.5738 0.4219 0.2728 0.7139 0.4972 0.2805 
0.4 1.2444 0.9033 0.5668 0.8267 0.6155 0.4092 0.9137 0.6428 0.3736 
0.5 1.6813 1. 23% 336 0.7990 1.2216 0.9159 0.6259 1.2229 0.8683 0.5178 
(pliant npenie fl gieecinorsinsanoninonin Pu Se TIERS TMT CETIO MBIT | : 
n=0 0.2213 0.1526 0.0839 0.2500 0.1700 0.0900 0 0 0 P 
0.1 0.2413 0.1675 0.0938 0. 2626 0.1791 0.0958 0.00898 0.00663 0.00428 
0.2 0.3030 0.2137 0.1244 0.3034 0.2089 0.1146 0.0388 0.0288 0.0189 
0.3 0.4164 0.2977 0.1793 0.3833 0.2683 0.15385 0.0994 0.0744 0.0499 
0.4 0.5977 0.4312 0.2673 0.5251 0.3745 0. 2256 0.2118 0. 1602 0.1110 
0.5 0.8787 0.6380 0.4036 0.7698 0. 5626 0.3619 0.4145 0.3190 0.2321 
} P13 Rin Pru . oF Ps — uM: 
»=0 0.4000 0.2667 0.1333 0 0 0 0. 1667 0. 1096 0.0524 a 
0.1 0.4180 0.2793 0.1405 0.01314 0.00943 0.00574 0.1732 0.1141 0.0548 | 
0.2 0.4766 0.3207 0.1649 0.0564 0.0406 0.0251 0.1947 0.1288 0.0629 | 
0.3 0.5917 0.4039 0.2154 0.1447 0.1054 0.0659 0.2368 0. 1582 0.0792 
0.4 0.7950 0.5518 0.3088 0.3067 0.2255 0.1456 0.3106 0.2103 0.1101 
0.5 1.1456 0.8126 0.4822 0.5976 0.4459 0.2992 0.4366 0.3012 0.1655 
Pi ~ 
“ = 0 0 0 0 SS. 
0.1 0.00481 0.00335 0.00190 
0.2 0.02055 0.01436 0.00829 ( 
0.3 0.0526 0.0373 0.0217 ( 
0.4 0.1112 0.0793 0.0478 
0.5 0.2155 0.1559 0.0965 ( 
do = (Pid + Foo + Foi) — + (11) Cp/oo = Pyb9 + P2b:A + P3619 + Pb:61 + (Psb2 — Nt 
Fa + Fé 4 Fb (12) Pa)d? ote (P75, <> P32)d9 +. (P62 wo Pa)dO, + 
1 = < 2 
eT ae - hie (Prd. — Pr2a)0o? + (Pisd, — Pya)000: + Eve 
. + + ‘ iw \A.2 5) 
by = y1(FrA + Fo + Foi) + Fit + Fura (18) (Pisde Pyea)61° (15) 
. ; ; ‘ Numerical values given in Table 3. che 
Numerical values for F; to F;; are given in Table 1. umerical values for P; to Py. are given in Table 3 
er 
; Longitudinal-Force Coefficient ciet 
Thrust Coefficient , iiss f 
The procedure used to evaluate the longitudinal- - 
The procedure used to evaluate the thrust coefficient force coefficient expression is identical to that used by fac 
expression 1s identical to that used by Wheatley. Sissingh,’ except that an additional term involving the I 
The simplified expression is coefficient of radial drag 6, was included. This term is tha 
cal 
C,/o0a = Ti + To + Toh (14) A(C,y/o0) = 2n76,/4 fo’? fo'cos* ¥(1 + tx) dy dx (16) v1 
Numerical values for 7 to 7; are given in Table 2. 6, is equal to twice the coefficient of skin friction for gg 
flow along the blades in the radial direction. The tice 
Power Coefficient final expression is cre; 
, bla 
The procedure used to evaluate the power coefficient Cy/oo = Hido + Heb; + Hb:0o + H45:0, + (Hobe + pov 
expression is identical to that used by Bailey®*® in Ha)? + (Hyd, + Hea) + (Hobe + Ha); + tior 
evaluating torque coefficient. The simplified expression (71152 + LTy2a)0o? + (isdg + Hygt)000, + F 
is (Ti552 + Hyea)0;* + Hyd, (17 tap 
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TABLE 4 
Numerical Coefficients for Equation for Longitudinal-Force Coefficient 
; 0 —0.4 —0.8 0 —0.4 —0.8 0 4 whe 
a aEn tata,’ | poerndni=<iee ieee cosedadpuagerulcedeaperenana 
n= 0 0 0 0 0 0 0 0 0 0 
0.1 0.0501 0.0368 0.0235 0.0147 0.0135 0.0120 0.0043 0.0088 0.0032 
0.2 0.1010 0.0743 0.0476 0.0305 0.0280 0.0253 0.0091 0. 0080 0.0069 
0.3 0.1534 0.1133 0.0732 0.0484 0.0445 0.0405 0.0151 0.0134 0.0116 
0.4 0.2080 0. 1544 0.1008 0.0695 0.0641 0.0586 0.0227 0.0204 0.0179 
0.5 0. 2656 0.1983 0.1309 0.0950 0.0877 0.0808 0.0330 0.0298 0.027 
ee ene a EE AERO | aes 
n=0 0 0 0 0 0 0 0 0 . 0 
0.1 0 0 0 —0.0559 —0.0339 —0.01389 0.1066 0 0763 0.0497 
0:2 0.0006 0. 0006 0. 0006 —0.1081 —0.0645 —0.0236 0.2153 0.1434 0.0996 
0.3 0.0021 0.0019 0.0017 —0. 1534 —0. 0872 —@Q.0241 0.3287 0.2352 0. 1505 
0.4 0.0050 0.0046 0.0042 —0.1875 —(). 0983 —(0.0098 0.4513 0.3213 0.2045 
0.5 0.0106 0. 0097 0.0090 —0. 2073 —0.0941 +0.0230 0. 5886 0.4197 0.2685 
—————_—_—_—_—_———H,—____—-~ ——————- -———H, -— NSN ————_—__—_—_—-H, ——----—-— — 
w=0 0 0 0 0 0 0 0 0 0 
0.1 —0.1126 —0.0808 —0.0502 0.2239 0.1701 0.1186 —0. 1036 —0.0755 —0.0475 
0.2 —0.2347 —0. 1683 —0.1041 0.4689 0.3568 0.2499 —0.2084 —0. 1504 —0.0924 
0.3 —0.3680 —0.2617 —0. 1587 0.7546 0.5750 0.4036 —0.3166 —0. 2252 —0. 1336 
0.4 —0.5179 —0 3639 —0.2145 1.1091 0.8458 0.5993 —0.4301 —0.2994 —0. 1698 
0.5 —0.6922 —0.4816 —0.2753 1.5722 . 2068 0.8724 —0.5530 —0.3765 —0.1997 
Hy——— ~ — —_——- HI, —_______“—~ ————_—__—_—_——— H,,_—_____—~ 
n=0 0 0 0 0 0 0 0 0 0 
0.1 0.1738 0.1298 0.0863 —0.0423 —0.0297 —0.0172 0.0891 0.0657 0.0425 
0.2 0.3593 0.2680 0.1767 —0.0892 —0.0629 —(0.0369 0.1938 0.14385 0.0943 
0.3 0.5714 0.4251 0.2798 —0.1474 —0.1045 —0.0625 0.3315 0.2484 0.1670 
0.4 0.8286 0.6143 0.40562 —0. 2293 —0.1650 —0.1032 0. 5307 0.4026 0.2802 
0.5 1.1557 0.8597 0.5739 —0.3544 —0.2617 —0.1754 0.8327 0.6435 0.4706 
——_—_————- 3 -—- SS Bias ———— — —Hy- a ai tn SS -Mis = 
p= 0 0 0 0 0 0 0 0 0 0 
0.1 —0.0700 —0.0489 —0.0278 0.1302 0.0933 0.0564 —0.0255 —0.0174 —0.0092 
0.2 —0.1481 —0.1037 —0.0594 0.2819 0.2031 0.1246 —0.0540 —0.0368 —0.0196 
0.3 —0.2442 —0.1720 —0.0998 0.4815 0.3510 0.2192 —0.0893 —0.0612 —0.0328 
0.4 —0.3747 —0.2660 —0. 1584 0.7668 0.5640 0.3646 —0.1362 —0.0934 —0.0511 
0.5 —0.5639 —0.4060 —0.2509 1.1941 0.8919 0.6001 —0.2022 —0.1396 —0.0776 
—-Hy.——_____—~ iy ———_~ 
p=0 0 0 0 0 0 6 
0.1 0.0475 0.0332 0.0188 0.0042 0.0034 0.0025 
0.2 0.1027 0.0719 0.0412 0.0170 0.0136 0.0102 
0.3 0.1749 0.1238 0.0719 0.0382 0.0306 0.0229 
0.4 0.2771 0.1974 0.1183 0.0679 0.0543 0.0407 
0.5 0.4281 0.3096 0.1914 0.1061 0.0849 0.0637 
Numerical values for H; to Hj; are given in Table 4. blades and used in the numerical coefficients. These 


Evaluation of Numerical Coefficients 


Inspection of the complete expressions for the rotor 
characteristics given in the Appendix will show that the 
terms that are to be converted into numerical coeffi- 
cients depend only on the advance ratio yu, the taper 
factor t, the blade mass coefficient y;, and the tip-loss 
factor B. 

Bailey has shown, and this analysis has corroborated, 
that the effect of large variations in ; upon the numeri- 
cal coefficients is small enough so that a mean value of 
¥1, substituted in the expressions for the numerical 
coefficients, will give satisfactory accuracy for all prac- 
tical blade designs. Other things being equal, an in 
crease in y, (corresponding to an increase in average 
blade coning) will cause a slight increase in the overall 
power required to maintain a given flight condi- 
tion. 

For this analysis, typical values for y; for the various 
taper ratios were estimated from data on existing rotor 


values were: 


—0.8 
14.41 


t 0 —0.4 
vy 5.00 7.52 


Similarly, a value of B = 0.97 was used throughout, 
since variations from this figure for normal rotor de- 
signs should be small. 

The numerical values were determined for a range of 
advance ratio » from 0 to 0.50 in 0.10 increments. 
Three values of the taper ratio, ¢ = 0, —0.4, and —0.8 
were substituted for all coefficients. 

For intermediate values of u, the coefficients may be 
plotted versus uw and the desired values read from the 
curves. For intermediate values of ¢, it has been found 
that linear interpolation is generally satisfactory. 


Derivation of Climb-Ratio Expression 


Consider the general case of a single rotor helicopter 
in steady, trimmed flight. A diagram of the forces 
acting is shown in Fig. 1. For convenience, the rotor 
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Fic. 1. Force diagram for helicopter in steady trimmed flight. 


center and aircraft c.g. are shown coincident, with the 
force vectors passing through the point of coincidence. 
In general, for an actual case, each of the forces shown 
would contribute a pitching moment about the aircraft 
c.g., the sum of which would be balanced by use of 
cyclic feathering control, deflection of a horizontal tail 
surface, or similar means. It is assumed that any 
additional external forces introduced for the purpose of 
trimming the helicopter are negligible. 

A diagram of the relative velocity components, to- 
gether with the angular relationships between these com- 
ponents and the helicopter flight attitude is shown in 
Fig. 2. The rate of climb is —V sin (a + 1). The 
forward velocity is V cos (a + +), which may be satis- 
factorily approximated by poR. 
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Fic. 2. Velocity diagram for helicopter in steady trimmed 
flight. 


The following assumptions are made: 

(i) T = T cos, since is a small angle. 

(2) Hsin c is negligible compared to 7, since H is 
a relatively small force. 

(3) D,sin (a + ¢) is negligible compared to 7, since 
whenever sin (a + ¢) is appreciable, D, is small, and 
vice versa. 

For steady flight, with the above assumptions, 


T=W | (18) 
In cases where an appreciable lift force is developed by 
nonrotor parts of the helicopter, this force should be 
subtracted from W in applying Eq. (18). 

A second condition for equilibrium in steady flight 
may be determined by stating that the sum of the force 
components along the flight path must equal zero: 

—T sina — Heosa+ Wsin (a + «) — D, =0 

(19) 


Solving for sin (a + +) and substituting Eq. (18), 
sin (a + «) = sin a + (H/T) cos a + (D,/T) 
(20) 

The climb ratio «x is defined as 
rateofclimb — —Vsin (a + «) 


= 21) 
rotor tip speed _ wR ' 





Substituting Eq. (20) in Eq. (21) and converting to coef- 
ficient form, 


k = —ptan a — wCy/Cr — w®Cp,/Cr cos’ a (22) 


This is the final form for the climb-ratio expression. In 
it, 
uw tan a = \ + C,/4B2(d\? + y?)”? (23) 


from Eq. (2) for a single-rotor-type helicopter, or its 
equivalent for other configurations. Similarly, minor 
modifications may have to be made in the foregoing 
derivation for configurations other than the single 
rotor. 


Construction of Climb-Ratio Charts 


A climb-ratio chart is a plot of climb ratio «x versus 
power coefficient Cp for a constant value of advance 
ratio uw. A set of such charts for values of » covering 
the speed range will express the complete performance 
of the helicopter, for any altitude, gross weight, and 
rotor speed. A typical example of a climb-ratio chart 
is shown in Fig. 3. As shown on the example chart, it 
will be found convenient to use blade pitch setting 
and thrust coefficient Cr as parameters. 

A suggested procedure for construction of such charts 
is: 

(1) Choose a desired range of values of 6) and Cr 
for the chart to cover the range of possible operating 
conditions for the subject helicopter. 

(2) Solve Eq. (14) for A, using combinations of the 
chosen values of 4 and Cr. 

(3) Using suitablé values of the physical constants 
do, 51, 52, 5,, 2, oo, and Cp, for the subject helicopter, 
form the required combinations for use in Eqs. (15) 
and (16) for C, and Cy. 

(4) Compute Cp, Cy, and «x for each of the chosen 
combinations of Cr and 6. The resulting values of 
Cp and x are plotted upon the climb-ratio chart. 
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Fic. 3. Example of a climb-ratio chart. 


Construction of Level-Flight Chart 


A level-flight chart is a plot of power coefficient 
required for level flight C, versus advance ratio yw. It 
expresses the complete level-flight performance of the 
helicopter, for any altitude, gross weight, and rotor 
speed. A typical example of such a chart is shown in 
Fig. 4. When a series of climb-ratio charts have been 
constructed for a range of values of the advance ratio 
u, the corresponding level flight chart may be con- 
structed simply by making a cross-plot of the climb- 
ratio chart values for x = 0. 


Blade Stall 


On the average helicopter there is a tendency for the 
tip region of the retreating blade to stall as high forward 
speeds are attained in powered flight. It is believed 
that the primary effects of such blade stall are excessive 


TYPICAL LEVEL-FLIGHT CHART 


| | 
—_—+—- t , ore a 





POWER COEFFICIENT Cp 





CE RATIO 


Fic. 4. Example of a level-flight chart. 


flapping of the blades and large variations around the 
cycle of the center of pressure location for those sections 
subject to stall, causing a sharp increase in rotor rough- 
ness and a tendency toward loss of control. The actual 
reduction in performance due to the energy loss in the 
stalled region should be small for a stalled region of 
moderate extent, since the resultant velocity at the 
blade will be relatively low. 

Even though excess power may be available, the ad-. 
vent of blade stall means that the helicopter is entering 
an undesirable flight regime that should be avoided. 

As an approximate criterion of the maximum blade 
section angle of attack experienced in a given powered 
flight condition, the angle of attack at the tip of the 
retreating blade may be calculated. It may be shown 
that this angle is 


= 0 + 0 + kg(ao + bi) + 
[A + a:)/(1 — w)] (24) 


Using this expression, limiting lines may be superim- 
posed on the climb-ratio charts and level-flight chart to 
indicate the boundary of the stalled blade region. A 
more complete treatment of blade stall is given by 
Bailey.® 


(a) y =3x/2, x=1.0 


DISCUSSION OF ANALYSIS 


Scope of the Complete Expressions for 

Rotor Characteristics 

The complete expressions for the rotor coefficients as 
derived for this paper include terms for the effect of an 
induced velocity distribution varying linearly from the . 
front to the rear of the rotor disc (X;), and the effect of 
use in the rotor design of means for varying the blade 
pitch angle linearly with the flapping angle (kg). These 
expressions are given in the Appendix in their unsimpli- 
fied form. (All A; terms have been dropped from the 
thrust, power, and longitudinal-force expressions, since 
examination of such terms has shown that they are 
negligible except in the expression for lateral flapping 
coefficient b;._ This is in accordance with the finding of 
Wheatley.*) 


Scope of the Simplified Expressions and Effects of 

Cyclic Feathering 

In simplifying the rotor characteristic expressions for 
use with numerical coefficients, it was found that un- 
desirable complexity resulted if the terms involving 
pitch change with flapping (k,) were retained; there- 
fore, such terms were dropped. It is believed that this 
restriction on the generality of the simplified expres- 
sions does not impair their validity for use in perform- 
ance calculations on helicopters incorporating pitch 
change with flapping, although in such cases the simpli- 
fied expressions for the flapping coefficients are, of 
course, not applicable. The reason for this belief is 
that the presence of pitch change with flapping means 
(within the assumptions of this analysis) that cyclic 
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feathering is produced by the flapping motion. This 
type of cyclic feathering is identical in effect to cyclic 
feathering produced by applied control. Lock? has 
shown that feathering is aerodynamically equivalent to 
flapping. Thus, the operation of a rotor subject to 
either or both types of feathering is equivalent to the 
operation of the same rotor as a pure flapping rotor (no 
pitch change with flapping or applied cyclic control) 
with its shaft axis inclined at some angle to the shaft 
axis of the actual rotor. If this proposition is correct, 
the only difference between a helicopter operating with 
feathering, and the same helicopter operating with a 
pure flapping rotor in the same trimmed flight condition, 
is a change in the fuselage and rotor shaft attitude with 
respect to the relative wind. In almost all cases the 
effect of such a change on the parasite drag, and thus 
on the performance, is negligible. 


Evaluation of Parameter K 


A recent theoretical investigation by the N.A.C.A. 
has provided a basis for the estimation of ‘the proper 
value for the parameter K in Eqs. (3) and (4). Asa 
very rough approximation, K may be taken as equal 
to unity for values of the advance ratio greater than 
» = 0.10. 


Effect of Blade Weight 


Bailey® has shown that for practical rotor blade de- 
signs the effect of the blade weight moment parameter 7 
on the rotor characteristics is negligible except for the 
flapping. Therefore, all terms containing 7 were 
dropped except in the expressions for the flapping coef- 
ficients. 


Offset of Flapping Hinge from Rotor Center 


The derivations have been made for a flapping hinge 
located at the axis of rotation. It is reasonable to as- 
sume that a displacement of the flapping hinge from the 
rotor axis a distance not greater than 5 per cent of the 
rotor radius will not impair the validity of the analy- 
sis. 


Use of Radial Blade Drag Coefficient 


In forward flight the rotor blades are, in general, 
operating in a condition of yawed flow. In order to 
account for the additional longitudinal force that must 
be present due to this effect, a contribution to H of the 
longitudinal component of the skin friction drag on the 
blades in the radial direction was included in the deriva- 
tion. The radial drag coefficient 6, is equal to twice the 
blade skin friction drag coefficient for flow in the radial 
direction. In most cases, it is probably sufficiently ac- 
curate to use 6, = do. 


Parasite Drag 


The parasite drag coefficient is intended to cover all 
parts of the helicopter except the blades. It is usually 


sufficiently exact to select a constant value of Cp, for 
use over the flight range. However, a variation with 
estimated fuselage angle of attack may be used if avail- 
able. 

Although a is usually a small angle, the cos* a term 
in the denominator of Eq. (22) should not be taken as 
unity. A simple chart may be prepared giving cos® « 
as a function of u tan a@ to simplify the determination of 
this term. 


Application for Hovering and Low Forward Speeds 


Results obtained by the method of this paper for 
hovering and low forward speeds are somewhat opti- 
mistic (see Fig. 5) in that the power required is under- 
estimated. This discrepancy is largely due to the basic 
assumption of a constant induced velocity distribution 
over the disc. This assumption introduces negligible 
error at medium and high speeds where induced power 
is a small fraction of total power, but at low speeds and 
in hovering, where induced power is of the order of 
70 per cent of total required power for the average heli- 
copter, the error becomes appreciable. 

The analysis should not be applied without suitable 
modification to the case of vertical or near-vertical 
descent at any but low rates of descent, since, in such 
cases, momentum theory conditions break down and 
the vortex ring flow state is entered. 


Autorotation 


The analysis may be applied to the autorotational 
(power-off) flight condition simply by extending the cal- 
culations into the range of low pitch settings. The 
curves on the climb-ratio and level-flight charts will 
then intersect the Cp = 0 line, thus covering the auto- 
rotational range. 


Application to Flight-Test Data Reduction 


The method of this paper may be used in the reduc- 
tion of helicopter performance-flight-test data. Test 
points may be reduced to coefficient form and compared 
with a set of charts calculated for the subject heli- 
copter. In this manner it should be possible to deter- 
mine consistent trends of differences between test 
and calculated performance, so that standard atmos- 
phere performance may then be predicted accu- 
rately. 


DISCUSSION OF FIGURES 


Figs. 3 and 4 


These figures have been discussed in the Analysis 
section. It may be suggested that for ease in deter- 
mination of coefficients, a simple auxiliary chart may 
be prepared for a specific helicopter which relates the 
nondimensional quantities used on the charts to their 
physical counterparts. 
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Fic. 5. Comparison of calculated and flight-test performance of 
a single-rotor helicopter. 














oO 


























































































































Y, 
L cv 
~ z 
ra S=0100-0250+.5000%| LY 
6 §-c110-0250+.50008 yh 
02 yee 
a 4, J 
wy LIA 
=! of 
5 LAT |NTEST DATA FOR 
sear = NACA 00l2 BLADE 
as os SPECIMEN, WOOD 
S L.E.,FABRIC COV- 
6 ERED, R.N.=172x 
ls = 
" 
ob 1, 





Yd a 6 i IZ 
ANGLE OF ATTACK &, DEG, 

Fic. 6. Comparison of analytical expressions for blade pro- 
file drag coefficient with wind-tunnel test results for a practical- 
construction rotor blade model. 


Figs. 5 and 6 


Fig. 5 presents a comparison between calculated and 
measured performance for a well-known helicopter. 
The flight-test points are taken from an excellent set 
of data collected by the N.A.C.A. during tests of the 
example helicopter, and published in a recent report. 
Values of parasite drag were available from wind-tunnel 
measurements. It may be noted that, except for hover- 
ing, agreement is good between the test points and the 
calculated curve for the higher assumed minimum pro- 
file drag value. 

Fig. 6 shows a comparison between.the analytical 
curve for blade profile drag coefficient used in the cal- 
culations for Fig. 5, and values found by the N.A.C.A. 
in wind-tunnel tests of a blade sample, of the construc- 
tion used for blades of the test helicopter. The higher 
analytical curve might represent the characteristics of 


a wood and fabric blade under service conditions, with 
roughened leading edge. 


Fig. 7 


Fig. 7 shows a comparison between climb perform- 
ance calculated for a typical small helicopter by the 
method of this paper and the performance as calculated 
by an approximate method.'® It may be noted that 
agreement is good in the low-speed range, but that the 
energy results are quite optimistic at high speeds. All 
possible factors in the calculations were kept the same. 
The blade profile-drag-coefficient variation for the cal- 
culations using the present method was chosen so as to 
equal the mean value of 6 = 0.0125 used for the energy 
method, at the mean C, of 0.66 as calculated by Wald’s 
expression. 


Fig. 8 


Fig. 8 shows the effect upon climb performance at a 
constant true air speed (approximately cruising speed) 
of the example helicopter of Fig. 7 due to major changes 


WALD'S ENERGY METHOD 
(REF. 10) 





0 





T » MPH. 


Fic. 7. Comparison of constant-power climb performance as 
calculated by this method and an energy method for a typical 
small helicopter. 


Data: W = 1,500 lbs. 565 = 0.0088 
R = 168 ft. 6, =—0.025 
f = 3.75 ft.? 6, = 0.590 
Cr/oo = 0.20 shaft hp. = 81 
B = 0.97 tip speed = 319 ft. per sec. 
a = 5.75 sea level, rectangular blades, ze1o 


twist 






» FUMIN. 
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RATE OF CLIMB 


ROTOR SHAFT HORSEPOWER 


Fic. 8. Effect of various changes on climb perfoimance of heli- 
copter of Fig. 7 at constant air speed of 87 m.p.h. 


@ f = 7.5 ft.2, d = 0.0088, wR = 319 ft. per sec. 
2) f = 7.5 ft.2, dd = 0.0088, wR = 425 ft. per sec. 

f = 75 ft.2, 8 = 0.0044, wR = 319 ft. per sec. 
 f = 3.75ft.2, 5 = 0.0088, wR = 319 ft. per sec. 
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in some of the helicopter parameters. The base curve 
(No. 1) corresponds to the performance of the heli- 
copter of Fig. 7 with its parasite drag area increased to 
7.5 sq.ft., a value that corresponds to present practice. 
Curve No. 4 illustrates the large benefits that may be 
achieved by halving the parasite drag. Such a reduc- 
tion is not unreasonable in view of the high parasite 
drag that is characteristic of most present helicopters. 
Curve No. 3, which might represent the use of the best 
available low drag’ section for the blades, indicates a 
considerably smaller, but by no means negligible, gain. 
The improvement shown by this curve is probably too 
great because of the sharp rise in profile drag coefficient 
at angles of attack beyond the low drag range that is 
characteristic of low drag sections. Curve No. 2 il- 
lustrates the effect of a considerable increase in rotor 
speed. The change in performance in the operating 
range is extremely small. This is contrary to the effect 
that would be produced in hovering, where the increase 
in rotor speed would increase considerably the required 
power. An increase in rotor speed is an effective 
means of reducing the tendency toward blade stall at 
high speeds. 





TRUE AIRSPEED, MPH. 
Comparison of additional shaft horsepower required 
for climb with rate of gain of potential energy in climb for 
helicopter of Fig. 7. 


Fic. 9. 


@ Power required for level flight. 

@) Additional power required for climb shown in Fig. 7. 

(3) Rate of gain of potential energy in climb (horsepower 
units). 


Fig. 9 


Fig. 9 shows an interesting comparison between the 
additional power required to climb and the rate of gain 
of potential energy in a rated-power climb as a function 
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of forward speed for the example helicopter of Fig, 7, 
It may be noted that the cost of climbing is less than the 
gain at hovering and low forward speeds, while the op. 
posite is true at higher forward speeds. The hovering 
and low speed effect is easily explainable on the basis of 
momentum and energy considerations. The effect at 
medium and high speeds is an indication that a lifting 
rotor becomes less efficient as it is required to contrib. 
ute more and more of a forward component of thrust 
as well as to continue to support the weight of the air- 
craft. It is for this reason that, other things being 
equal, an autogiro or helicogiro (in which all or part of 
the engine power is delivered to a conventional propel- 
ler) would probably have a slightly higher top speed 
than the equivalent helicopter,-in spite of the power loss 
due to a propeller efficiency less than unity. 
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Appendix 


Flapping Coefficients 


The complete forms of the three simultaneous equations that determine the flapping coefficients dp, a1, and b; are 


0 = Al!/sB* + 0.07953 + t(4/,B* + 1/seu*)] + (00 + Redo) ['/sB4 + */,B2u? ~ */sou* + £(1/sB5 + 1/6B%u*)] + 


0:[1/sBS + */eBy? + t(1/eB* + 1/sB4u*)] — kabi[*/sB%u + t('/<B'u)] — (1/n)(@o + 17) (25) 


0 = X[/eBYy — au? + t(/sB%u — 0.053u')] + (0 + kpdo)l?/sB%x + 0.053u4 + t(1/2B%)] + 61["/2B + t(2/sBu)] + 
a,[—1/,BS + */sB%u? + t(—1/sB5 + 1/iBu2)] — kpbi[!/B4 + 2/sB%u? + £(1/sBS + 1/,B%u*)] (28) 
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0 = M/eBY + £(/sB)] + aol4/eBYe + 0.035u* + £(1/4B%)] + haar ["/4B* + 4/eB%? + t(2/sB + 1/rpBYu?)] — 
bs [1/.B* + 1/eBu? + t/.BS + 1/BYy2)] (27) 


Thrust Coefficient 


The complete expression for the thrust coefficient is: 


- = A{1/2Bt + 1/au? + [!/sB* + (2/9x)u*]} + (00 + kpao) [*/sB* + */2Bu* — (4/9m)u* + t(*/.B* + 1/,B%u* — */su*)] + 
¢ 
6:[1/.B4 +1/,B%p? — */sou* + t(1/sB5 + 1/gB%u2)] + a1(*/qu*) — bobs {*/2B%e + */au* + t["/sB%u — (4/453)u*)} (28) 


Power Coefficient 


The complete expression for the power coefficient is: 


C. 
- = 5o[*/4 + 47am? — '/sou* + t(*/5 + */eu?)) + 5 {A((*/s + */4t) + (0 + kpac)[*/e + au? + tC°/s + ou®)) + [1/5 + /eu? + 
a 


0 
t(/4 + */su)] — Rpdi(*/su + */aut)} + 52(A2{ 4/2 — */qu? + £[*/s — (4/9e)u*]} + Ao + pao) [2/2 + (4/9x) ut + (2/2 + */rou*)] + 
NA(1/2 + /ieu* + 3/st) + (00 + Rpao)?[*/e + 1/au? — */so4 + t(1/5 + 1/eu*)] + O1(00 + pao) [2/5 + 2/am? + t(1/2 + */qu?)) + 
O:2[?/6 + 1/su? + t(2/r + 1/10m?)) + air(?/eu — 2/qu® + 1/sut) + kpdid(—'/ou + */su® — 4/gut) + kpbi(Oo + kpao)(—*2/su — */out) + 
hpbi0i(—*/2 — */sut) + ao2(2/4u? — */ieu* + 1/eu%t) + ao(kpar — b1)(4/su + */aut) + a1?[*/g + 2/16? + t(4/10 + 1/eu*)] + b1?[*/8 + 
V/rgu? + t(*/1o + */eau®)] + Rpaidi(—*/qu? — */ou*t) + kp2(ai? + by) [*/s + */su? + t(*/10 + 2/eu*)]) — a(a?*{'/,B? — Vy? + 
t[!/sB* — (4/9x)u}} + (00 + pao) [1/s3B3 + (2/9) u* + £(2/4B4 + 1/sou*)] + Oi [1/4B* + 3/sou4 + £(2/sB5)] + ai d['/2B%u — 8/eu® + 
t(#/;B%u)] + ked[—*/eB® + */16m? + t(—1/sB*)] + ao? [*/4B2u? — 1/rgu* + t(1/6B3u*)] + ao(?/okprr — 0,)['/sB3 6 + t(1/.Btu)) + 
a;*[1/gB4 + 2/:6B2yu? + t(2/10B® + */gB%u2)] + b:2[*/gB4 + */1gB2u? + £(4/10BS + */o4B%u?) | + kpaydi[—'/sB%u? + t(—'/12B%u2)]) (29) 


Longitudinal-Force Coefficient 


The complete expression for the longitudinal-force coefficient is: 


= bof "au + '/au® + t['/am + (4/45x)u8]} + 8,[(4/3)u2(1 + /at)) + &{AC/am + Y/aut) + Oo + hpa0)(*/am + */aut) + 0:(*/su+ 
0 


1/sut) + ai[—*/6 + */eu? + t(—"/p + */iou®)] — kpbil*/e + */eu? + t(2/e + 3/rou?) ]} + 52(A7{ */ou + t[(2/3)u2]} + AO + keao)[ u— 
(4/3x)u? + t(*/eu — "/eu*)] + if */ou — 1/eu® + t[*/e — (8/45m)ut]} + (00 + kpao)*{*/en + */gu® + t[*/au + (4/450) u*]} + 
(00 + p20): [2/3 + (8/45m) u* + */out] + 0:2(4/au + */sut) + bear(A{ —1/2 + 5/gu® + t[—1/3 + (44/45) u*]} + (Oo + hea) ([—1/s + 
Way? — (28/45) u® + t(—2/4 + */eu?)] + Or[—*/e + 1/su® + H(—*/s + /i2m?)]) — Sekpi(A{*/2 + 8/eu? + t[*/s + (16/45) us!} + 
(00 + kBao) ["/s + 3/su? — (832/45m)y® + t(2/4 + 3/au)] + O:[*/4 + 3/eu? — (5/96) u* + t(4/s + */qu) ]) + 50(ao?{ */qu? +2 [(2/154)u*]} + 
ao(kpa, — 01) [*/au? — (4/154) uy? + 1/gu%] + 2/4[ar2(bp? — 1) + b,2(3kp? + 1)](*/on + 1/sut)) — a(A®["/ou + (2/3) u%] + Ao + 
kao) (!/2Bu — (2/3) u* + t(2/.B%e — */rgu®)] + Oi {1/4B2n — 8/1? + t['/:B% x — (4/45) u*}} + ad{ —3/,B? + %/seu? + t(—1/2B* + 
(14/159) u*]} + a1(00 + kaa) (—1/sB*— (10/45) u3—1/,Bt] + a10i(—1/4B* — 1/,Bt) — kabid{*/4B? + */rgu? + t[*/B* + (8/452) u*]} + 


ast{ —1/.BYu + */rgu + t[—1/eB%u + (4/45x)u*]} + ao(bi — haar) |"/eB* + (2/45m)u* + ¥/sBY] + ax(Robi — a1)(1/B%u + 1/eBut)) 
(30) 





Letter to the Editor 


pends upon the bending rigidities of the plate in the longitudinal 
and transverse directions as well as its torsional rigidity at the 
moment of buckling: 


Ger. = KE rea. (t/b)? (3) 


Dear Sir: 

The writer fully agrees with George Gerard’s statement in his 
paper entitled “‘Secant Modulus Method for Determining Plate 
Instability Above the Proportional Limit,” published in the Janu- 





ary, 1946, issue of the JouURNAL OF THE AERONAUTICAL SCIENCES 
that “the critical stress and critical strain are implicitly related 
by the stress-strain curve of the material, i.e., the secant modulus 
relates the critical stress and critical strain.’’ When the correct 
critical strain €,r. in plate buckling is known the critical stress ocr. 
follows immediately by multiplying ¢.. by the secant modulus 
Exec. There is no foundation, however, for the belief that the 
formula 

Ger. = KE(t/b)? (1) 


which obviously gives an incorrect value for the critical stress 
beyond the elastic range when Young’s modulus—valid only in 
the elastic range—is substituted into it, will predict the critical 
strain correctly above the elastic limit if written in the form 


éer. = K(t/b)? (2) 


In the inelastic range the value of E in the right-hand member 
of Eq. (1) must be replaced by a reduced value Erea. which de- 


The right-hand. member of Eq. (2) can be obtained from the 
right-hand member of Eq. (3) by dividing the latter by Epea., 
while the left-hand member of Eq. (2) can be obtained from the 
left-hand member of Eq. (3) by dividing the latter by Eseo. Ob- 
viously an equation loses its validity if its two members are di- 
vided by two different quantities. 

If the author’s reasoning were valid fox plates it should hold 
for columns also. Hence the buckling strain of a column in the 
inelastic range should be given by 


€or. = 4*/(L'/p)? (4) 
instead of by 
fer. = [4*/(L'/p)*] (Ered. /Enco-) (5) 
which formula follows directly from the von K4rm4n equation 
Cer. = WE rea. /(L'/p)? (6) 


(Continued on page 380) 








Analysis of Relay Servomechanisms 


HERBERT K. WEISS* 
Aberdeen Proving Ground 


SUMMARY 


Methods are developed for the transient and steady-state 
analysis of relay-type servomechanisms, with intermittently ap- 
plied torque, inertia, viscous friction, dry friction, finite spacing 
of the controlling contacts, constant delay time in torque change, 
and stabilization by ‘‘flip-flop,’”’ ‘‘front-lash,’’ and anticipatory 
circuits. 

Solutions are obtained through representation of the motion 
of the servomechanism in the position-velocity or ‘‘phase’’ plane. 
In this plane, the discontinuous character of the driving torques, 
dry friction, delay time, finite contact spacing, and certain types 
of stabilizers, defines boundaries that are frequently straight lines. 
The boundaries divide the phase plane into a number of regions, 
within each of which the motion of the servomechanism is de- 
terminable by elementary means. Theorems and semigraphical 
constructions are developed for the rapid representation of the 
transient and steady-state motion in the phase plane, and simple 
constructions are demonstrated for converting the phase-plane 
plots to conventional time plots. 

Guided by phase-plane representations, exact mathematical 
solutions are obtained for a number of cases of interest. A 
method is developed for testing the physical significance of the 
mathematically derived steady-state solutions. Finally, the 
techniques developed are applied to the comparative examination 
of several methods for suppressing ‘‘hunting’’ of the servo. The 
analysis indicates an anticipatory signal to be the most promising 
for control of torque application in a high performance relay 
servomechanism. 


INTRODUCTION 


I ONE OF THE EARLY PAPERS on the quantitative de- 
sign of servomechanisms, Hazen! discussed several 
broad classifications of servomechanisms, including the 
relay type and the continuous type. Because the latter 
type is amenable to analysis by familiar mathematical 
tools, ?~* its discussion has predominated in the litera- 
ture. Nevertheless, commercial and military applica- 
tion of the relay- or clutch-type servomechanism has 
proceeded,’~'° sometimes in situations where the relay 
servo matched the continuous type in performance and 
surpassed it in simplicity and light weight. Fry’ has 
given an excellent qualitative description of modern 
relay servos, which he calls “‘stepping follow-ups,’’ in- 
cluding a discussion of stabilizing methods. 

Operation of a relay servo is characterized by abrupt 
changes in the values of the fundamental parameters, 
primarily driving torque. Between these abrupt 
changes, the value of the parameter is substantially 
constant. In the case of torque, the discontinuities 
are introduced by control through relays, clutches, 
contacts, and other intermittently operating devices. 
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Two examples of relay-type servomechanisms de- 
scribed in the literature’—® are indicated in schematic 
form in Figs. la and 1b. One is a split-field series d.c. 
motor with the fields controlled by a relay. In the 
commercial version,’ the relay is, in turn, controlled 
by a vacuum tube amplifier. The second type is a 
clutch-type servo, with the direction of torque applica- 
tion controlled by solenoid-operated clutches. Also 
usable in this fashion are reversible shaded-pole induc- 
tion motors, reversible single-phase capacitor-type in- 
duction motors and, in fact, any power source with 
reversible output which can be controlled by a two 
position switch. It is this flexibility with regard to 
torque source which makes the relay-type servo an in- 
teresting possibility for many applications. 


METHOD OF ANALYSIS 


The fact that the nonconstant parameters related to 


the relay servo can, in many cases, be considered to be 
functions solely of servo output angle and velocity, or 
error angle and velocity, suggests transformation of 
the analysis to the position-velocity or phase plane. 
This plane has been.used successfully by G. R. Stibitz 
in an unpublished analysis by statistical methods of a 
pulse-controlled servomechanism, by DeJuhasz" for 
the analysis of nonlinear processes and instruments, 
and by Minorsky” for the study of circuits parametri- 
cally excited. 
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Fics. la and 1b. Schematic drawings of two typical relay-type 
servomechanisms. 
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ANALYSIS OF RELAY 


For an initial exploration of methods of solution in 
the phase plane, a simple, idealized relay servomecha- 
nism, identical to that of Hazen’spaper,'isassumed. Its 
elements, shown in Fig. 2, consist of a torque source, 
inertia, viscous damping, dry friction, and controlling 
contacts that have a finite inactive zone. The torque 
is considered to be plus or minus a constant value, or 
zero, depending upon the position of the controlling 
contact. Time lag is assumed to be representable as a 
constant delay time between the contact operation and 
the resulting torque change. It is further assumed that 
delay time in creating torque is equal to delay time in 
removing it. 
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Fic. 2. Elements of the idealized relay servo. 


The elements of the relay servo of Fig. 2 are identi- 
fied by symbols in the following tabulation. The num- 
ber of parameters is then reduced, and the generality of 
the processes improved, by introducing dimensionless 
parameters, after the fashion of Draper,'* who has em- 
phasized the value of this step in the analysis of instru- 
ments. 


SYMBOLS 


= moment of inertia of the servo, referred to its output 
shaft 
= viscous damping coefficient, referred to the output shaft 


c 

@, = angular position of the servo output shaft 

6; = angular position of the servo input shaft 

@. = angular inactive zone between the two contacts 

T = servo driving torque, referred to the output shaft 

F = dry friction torque, referred to the output shaft 

tg = constant time delay in application or removal of torque 
T after contact closes or opens 

t = time 


The following parameters are derived from those defined 


above. 
= [/c = “time constant”’ of the servo (1) 
6, = T/c = maximum servo output rate with one con- 
tact closed (runaway velocity of output) (2) 
5 = 60, — 0; = servo error (3) 
w = 6; when input rate is assumed to be constant (4) 
s = t/r = dimensionless measure of elapsed time (5) 
Sa = ta/r = dimensionless measure of time delay (6) 
o = 6/6,r = dimensionless measure of magnitude of 
error (7) 
a, = 6./6;r = dimensionless measure of contact spacing (8) 
ry = |F/T| = absolute value of ratio of dry friction 


torque to servo driving torque (9) 
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tw = |w/6,| = absolute value of ratio of constant input 
velocity to maximum servo velocity (10 


DIFFERENTIAL EQUATIONS OF MOTION 


When the inactive zone between the two controlling 
contacts is of negligible width, the sign of the applied 
torque depends upon the sign of the error #, — 4,. 
This dependence will be symbolized by the expression 


Torque = — 7 sgn (6, — 0;) = —TZ sgn (6) (11) 


The dry friction torque is expressed as 
Friction torque = —F sgn (6n) (12) 


For finite width of the inactive zone between con- 
trolling contacts, it is convenient to locate the zero 
value of 6 at the edge of the contact that drives the 
servo ina positive direction. Replacing the dimen- 
sional measure of error 6 by the dimensionless measure 
o allows the driving torque to be expressed as a function 
ot o 


Torque = —Tf(c) (13) 


where f(a) is —1.0 for ¢ < 0, 0 for 0 < o < a, and 1.0 
ford > a. , 

The differential equation of motion for the relay 
servo is 


16» + Cm = —Tf(c) —F sgn (6) (14) 
Introducing Eqs. (1), (2), (3), (4), and (9) 
78 +6 = —6,f(c) — w — rpb,sgn (6 +.) (15) 
Now differentiating Eq. (7) with respect to s 
(dé ‘dt )(dt/ds) 


bs = 6 /6, (16) 


dc/ds = § = 
Similarly 


da/ds = ¢ = (rg) /6, (17) 


By introducing Eqs. (16) and (17) to Eq. (15), the 
measure of time is changed from the dimensional unit ¢ 
to the dimensionless unit s. When this is done, and 
Eqs. (8) and (9) employed, the differential Eq. (15) 
becomes 


o+o= —If(o) —7r, — resgn(a+r,) (18) 


BOUNDARIES IN THE PHASE PLANE 


Since the functions f(¢) and sgn (¢ + 7,), which pre- 
vent a general solution of Eq. (18) from being obtained, 
depend only upon o and 4a, they define contours in the 
phase plane (position-velocity plane). These contours 
(straight lines) divide the plane into six regions. In 
each region Eq. (18) can be solved by elementary 
methods. Eq. (18) can be written in generalized form 
for any of the six regions 


coto=k (19) 


where & is a constant in that region, and is the sum of 
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Fic. 3. Regions of constant & in the phase plane. 
the terms to the right of the equal sign in Eq. (18). 
The form of Eq. (18) in each region and the shape of the 
regions are shown in Fig. 3. 


SOLUTION OF GENERALIZED EQUATION 


It is convenient to solve Eq. (19) for the initial values 


s=0,¢=0,¢=0 (20) 

whereupon 
é/k =1—e~“ (21) 
(o — oo)/k = s — (1 — e~*) (22) 


Following the nomenclature of Minorsky,'* a curve 
in the phase plane representing velocity-position points 
occupied successively by the servo will be called a tra- 
jectory of the servo in the phase plane. 

Eliminating s from Eqs. (21) and (22), 


oo/k = (o/k) + (a/R) + In [1 — (@/k)} = (28) 


From Eq. (23) a sheaf of trajectories can be con- 
structed for each region, describing the passage of the 
servo through the region, from any initial combination 
of position and velocity. Fig. 4 shows such a sheaf for 
k = 1.00. Note that the trajectories are identical, ex- 
cept for a horizontal shift along the o axis. This fact 
is also indicated by Eq. (23) and will facilitate the 
graphical solution of problems, since only one trajec- 
tory need be computed for each region. 


THEOREMS 


Two relationships derived from Eqs. (21), (22), and 
(23) will be of sufficient later use to be designated as 
theorems. They are: 

(1) Any two points (1) and (2), lying on the same 
trajectory are related by the expression 


























Fic. 4. Trajectory sheaf for particular value of 


(o1/k) + (o:/k) + In{l — (0:/k)] = (02/k) + (62/k) 
+ In[l — (62/k)] (24) 


derived from Eq. (23). 

(2) If two points (1) and (2), lying on the same 
trajectory, are separated by the time interval s. — 5, 
where 52 is the time corresponding to the second point 
and s; is the time corresponding to the first point, the 
rates at the two points are related by 


[1 — (G2/k)] = [1 - (o:/k)Je—“*-™ — (25) 


derived from Eq. (21). 
For the case of k = 0, Eqs. (24) and (25) reduce to 
0} + o1 = 02 + G2 (26) 
oe “ie oe— (27) 
and the trajectory in the phase plane becomes a straight 
line of slope —1. 


TIME MEASUREMENT FROM THE TRAJECTORY 


Since much of the following development will be con- 
cerned with graphical constructions in the phase plane, 
it will be useful to have a simple method for converting 
the deduced trajectories to more conventional plots of 
o and o against s. From Eqs. (21) and (22) the elapsed 
time s from ao to any point (1) on the trajectory is de- 
termined to be 


ksy = 01 + a1 — a (28) 


ks, is determined by a simple graphical construction 
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ANALYSIS OF RELAY 


in the phase plane (see Fig. 5a). The intercept on the 
g axis between o» and the intersection of a 45° line from 
point (1) is k times the elapsed time from oo to point 
(1). The intercept on the o axis between a 45° line 
from point (1) and a 45° line from any other point on 
the same trajectory (2) is k times the time difference 


Ss. -- S, between these two points. Since k changes 
from region to region, an additional construction to 
eliminate k and obtain s, — s, directly is desirable. 


The axes ¢ = kand o = 1 — kwill always be available. 
Fig. 5b shows how, by the addition of one line, the time 
difference sz — s, is determined directly as an intercept 
on the 1 — & axis. 
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Fics. 5a and 5b. Graphical determination of time differences 








The graphical construction fails when k = 0, but the 
time difference can then be obtained from Eqs. (26) and 
(27), noting that ’ 


tan B = ee) —l= (d2 = a1) /o2 (29) 


where 6 is an angle easily measured in the phase plane. 


CONTOUR OF CONSTANT TIME RELAY 


A third theorem will now be stated, defining a general 
relationship existing between events on all trajectories 
in the same region for a particular set of initial events. 

(3) Ina particular region, if the locus of points cor- 
responding to a given event happening on each tra- 
jectory is a straight line, described by 


oi/k = a + b(6;/k) (30) 
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then the locus of events on each trajectory, following 
the initial event by a constant time delay s, is also a 
straight line, and is described by 


o2/k = A + B(a2/k) (31) 

where 
A =a+t sq — (1 + db)(e* — 1) (32) 
B= e4(1 +b) -—1 (33). 


This theorem is readily derived from Eqs. (24), (25), 
and (30). 


TRAJECTORY WEB FOR SERVO WITH CONSTANT 
VELocITY INPUT 


An assembly of trajectories traversing the various 
regions of the phase plane will be called a trajectory 
web. Such a web has been constructed for a relay 
servo with inertia, viscous damping, contact spacing, 
and constant velocity input. It is shown in Fig. 6. 
Arrows indicate the direction of motion of the servo 
along the trajectories in each region. Regardless of 
the initial ¢, ¢ point that may be chosen, tracing of the 
trajectories from region to region eventually leads to 
the origin, confirming Hazen’s deduction! that the 
steady-state error of a servo of the type assumed de- 
pends only upon the contact spacing, for constant veloc- 
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Fic. 6. Trajectory web for relay servo with inertia, viscous 
damping, finite contact spacing, and constant velocity input. 
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EXAMPLE OF GRAPHICAL SOLUTION WITH TIME LAG 


Graphical determination of the transient and steady- 
state motion of a relay servo with inertia, viscous damp- 
ing, time lag, and constant velocity input has been car- 
ried out to demonstrate the principles developed thus 
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Fic. 7. Phase plane solution for the transient and steady- 
state motion of a relay servo with inertia, viscous damping, 
constant delay time, and constant velocity input. 


far. Hazen! has shown that a servo with these charac- 
teristics hunts in its steady state. Two initial points 
have been chosen, the first with the initial error ex- 
tremely large, and the second with the initial error ex- 
tremely small. The trajectory web, Fig. 7, shows how 
the servo motion, in both cases, converges on a continu- 
ous hunting state. Although successive full cycles of 
the oscillation in the steady state are identical, the half- 
cycles are, unlike those of a continuous linear servo, 
unsymmetrical about zero error. 

Time lag shifts the torque interchange contour in 
Fig. 7 from the o axis to the inclined axes. These axes 
are determined from Eq. (31) by setting a and } equal 
to zero. 

Magnitude and velocity of error at peak amplitudes 
can be read directly from Fig. 7. Elapsed time between 
these peaks can be determined by the graphical con- 
struction given earlier. Variation of magnitude and 
velocity of error with time can then be plotted, and this 
has been done in Fig. 8. Note how different these 
curves are from the familiar sine waves of continuous, 
linear servomechanisms. 


TimE LAG AND CONTACT SPACING 


When contact spacing, time lag, and constant veloc- 
ity input are simultaneously considered, it is found nec- 
essary to standardize on a subscript system, Each 
contour is numbered. Points lying on a contour bear 
the same number as the contour, and the oa, o coordi- 
nates of the points bear the subscript of the contour 
that they define. Fig. 9 shows a servo path circling 
the origin. At (1) the positive driving torque is ap- 
plied. The servo crosses the o axis at (0), and at (2) 
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Fic. 8. Variation of error and error velocity with time for a 
relay servo with inertia, viscous damping, constant delay time, 
and constant velocity input. 
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Fics. 9a and 9b. 
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the positive driving contact opens. At (3) the negative 
contact closes. At (4), which is sg later in time than 
(2), the positive driving torque is removed, and at (5), 
which is sg later in time than (3), the negative torque is 
applied. The sequence is repeated as the servo passes 
points (6) to (1). 

Fig. 9a is drawn for the case 


Sa > Se (34) 


where s, is the time required for the servo to cross the 
open contact interval o,. A similar figure, but for 


$a < 3 (35) 


is presented as Fig. 9b. Note that since points (2) and 
(3) are s, apart in time, so likewise are points (4) and 
(5). This fact is made use of in solving for the (5) 
contour. 


GRAPHICAL SOLUTION 


When sg < s,, contours (5) and (1) are obtained di- 
rectly from Eq. (31) with appropriate values of a, }, 
and k inserted. When sg > s,, Eq. (31) is no longer 
valid, as the servo path from (3) to (5) passes through 
two k regions. A graphical constfuction in the phase 
plane is the simplest way of determining the (5) con- 
tour. Applying Eq. (25) to points (2), (3), (4), and 
(5), and eliminating s,, 
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The corresponding expression for points (7) to (1) is 
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T¥OVee (37) 
Note that the friction term rp has been dropped in the 
above expressions. The efforts of the servo designer 
will almost always be directed toward elimination of 
dry friction in a high-performance relay servo, and the 
comparatively small importance of the rp term is not 
considered sufficient to warrant carrying it through 
further developments. In any practical case where 
the dry friction of the load is important, the servo 
analyst will find no difficulty in restoring the rp factor 
to the analysis. 

The form of Eqs. (36) and (37), in which the only 
unknowns in the graphical analysis are o, and a, sug- 
gests a solution by similar triangles. The necessary 
parameters are found to be located conveniently on the 
plot in the phase plane. The similar triangles, ab- 
stracted from the phase plane for clarity, are shown in 
Fig. 10a. The construction, as actually carried out in 
the phase plane, is shown in Fig. 10b. 

As an example of the use of Eqs. (36) and (87), a 
trajectory web has been constructed for a relay servo 
with inertia, viscous damping, time lag, contact spac- 
ing, and constant velocity input for values of the para- 
meters such that the steady-state oscillation takes 
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Fic. 10. Graphical solution for determining contour of applica- 
tion of opposing torque. 


place with sg > s,. The web is shown in Fig. 11. The 
values of sz and o, used have been chosen for clearness 
of presentation and do not necessarily indicate typical 
values for any real servo. 


MATHEMATICAL SOLUTIONS 


The graphical methods of analysis developed so far 
furnish a rapid and simple means for the study of par- 
ticular relay servos. At the expense of considerably 
greater computational time, surveys of the steady-state 
motion, as affected by parameter variation, can also be 
accomplished mathematically. The trajectory plots 
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Fic. 11. Phase plane solution for the transient and steady- 
state motion of a relay servo with inertia, viscous damping, 
constant delay time, finite contact spacing, and constant velocity 
input. 
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in the phase plane, however, are still found to be invalu- 
able aids in guiding the mathematical representation to 
a solution. 

Mathematical determination of the steady-state am- 
plitude of hunting of a relay servo with inertia, viscous 
damping, time lag, and contact spacing will now be 
made. Input will be assumed stationary. Solution 
is accomplished by assuming first that the servo does 
hunt in the steady state and solving for the amplitude of 
oscillation. Physical existence of this oscillation is then 
verified by ascertaining whether or not the servo will 
return to the conditions of the solution if differentially 
disturbed from them, An alternate, and quicker, veri- 
fication can be accomplished by graphical construction 
of the case under consideration in the phase plane. The 
graphically represented servo, behaving as a true ana- 
logue, will not converge to any steady-state oscillation 
that has only mathematical existence. 

The case of sz > 5s, will first be examined. Since in- 
put is assumed stationary, only one-half cycle of os- 
cillation need be considered. Refer to Fig. 9a. Points 
(1) and (4) are related by Eq. (24). Points (2) and 
(4) are related by Eq. (25). Points (4) and (5) are 
related by Eqs. (26) and (27). o4 and o are related to 
each other by Eq. (31) with a and 6 zero, and the same 
expression relates og to 9. Assumption that a symme- 
trical steady-state oscillation exists provides the con- 
tour-closing relations 


= “th + oe (38) 
os = —o1 (39) 


Employing only the designated equations, and elim- 
inating all o and ¢ among them except 02, the relation is 
obtained that 


e*4[(1 — o2)e*~*] — e4(1 + e~) + e7 (1 — o2) = 0 
(40) 

The maximum velocity reached in the steady-state 
oscillation is not g2 but os. Algebraic shuffling of the 
equations indicates, however, that a2 is more easily cal- 
culated than oy. o« is readily determined from o2 by 
Eq. (25). 

When contact spacing is zero, Eq. (40) reduces to 
Hazen’s expression’ for the g2 velocity of steady-state 
oscillation of a relay servo with time lag and zero con- 
tact spacing. 

In order to solve Eq. (40) numerically, it is necessary 
to eliminate s,. Applying Eqs. (24) and (25) across 


Sey 


(l—w)(l—e") +o,—s,=0 (41) 


From Eqs. (40) and (41) the effect of contact spacing 
o, and time lag sz, on velocity 2 can be determined, and 
thus the effect of these parameters on maximum veloc- 
ity in hunting o«. The solution, however, is good only 
for sg > 5S,. 
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When sg < s,, examination of Fig. 9b reveals that 
Eq. (40) is still valid but Eq. (41) must be replaced by 
the following expression 

Sat oz—o,—e-“(e%4#—1+62)=0 (42) 
which may be combined with Eq. (40) to give 


= (sa — % + a2) 
eSa—eeFs) (982/7 oe a2)] ms 





e*¢ (43) 
Eq. (42) was obtained by applying Eqs. (24), (25), 
(26), and (27) across the o, interval. 


RESULTS OF COMPUTATION 


Fig. 12 shows the effect on o2 and go of variations in 
contact spacing ¢, and time lag sz. The boundary s, = 
$s, between the two regions for which solutions were ob- 
tained is shown. The amplitude of swing outside the 
contact inactive interval is also given. This is com- 
puted from Eq. (23) to be 


2 = o2 + In (1 — oz) (44) 


The full double amplitude of oscillation o, including 
the contact spacing is 


a, = 20, + o (45) 


Note in Fig. 12 that for some values of sz and o, there 
appeat to be two solutions for amplitude of oscillation 
o,. Attempts to observe the smaller amplitude in the 
phase plane by graphical methods fail to reveal this 
solution as a steady state reached after a transient dis- 
turbance. Hence, a mathematical determination of the 
stability of the solution is desirable. 


STABILITY OF SOLUTION 


The solution will be called stable if it represents a 
continuous oscillation to which the servo returns after 
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Fic. 12. Amplitude of steady-state oscillation of a relay servo 
with inertia, viscous damping, constant delay time, finite contact 
spacing, and stationary input. 








ils that 
uced by 


(42) 


(43) 


), (25), 


ions in 
Y Sa = 
pre ob- 
de the 
$ com- 


(44) 
luding 


(45) 


, there 
lation 
in the 
il this 
it dis- 
of the 


nts a 
after 


e 


aL 


ITY 
No 


ervo 
itact 





ANALYSIS OF RELAY SERVOMECHANISMS 371 








OF 
5), 
tb 
2 


tet 
Pg 











<< 


Fic. 13. Trajectories for determining stability of a solution. 


small disturbances. For example, suppose that o,, 
o, are the values given by the solution for the point at 
which the servo trajectory, in its steady-state motion, 
crosses contour (4). Assume, however, that the servo 
actually crosses contour (4) with a velocity o4, + Ado,, 
where Ag, is a small quantity. This is shown in Fig. 
13, where contour (A) is that defined by the solution, 
and contour (B) is differentially displaced from contour 
(A). From (4), the servo oscillates, its trajectory in 
the phase plane circling the origin, and again crossing 
the (4) contour with a velocity os, + (Ags)o. (Ads)e 


is the deviation from 6,4, after one complete cycle. 


Now, if (Agy)e < Ady, it is evident that each successive 
cycle will bring the displaced trajectory closer to the 
trajectory given by the mathematical solution, and the 
latter correspondingly represents a continuous oscilla- 
tion, in which the servo will persist indefinitely. More- 
over, the amplitude and period are perfectly definite 
quantities, to which the servo will return, if differen- 
tially displaced by external influences. Here again, the 
relay servo differs fundamentally from the linear con- 
tinuous type, for, if the latter is adjusted to hunt in- 
definitely, the amplitude of steady-state oscillation .de- 
pends upon initial conditions, and can be changed at 
will, 

One finds in practice, of course, that no continu- 
ous servo is perfectly linear, and the variation in param- 
eters generally assumed constant can cause a continu- 
ous servo to have a stable amplitude of hunting, similar 
to that of the relay servo. 

When, for the relay servo, (Ags). > Ady, each cycle 
carries the servo farther from contour (A), after an 
initial differential displacement, and the solution given 
by the mathematics represents a state that cannot be 
maintained by a physical servo. 


DETERMINATION OF SOLUTION BOUNDARIES 


Examination of Fig. 12 indicates that the sg < s, 
region is the most promising one in which to study the 
stability of the results depicted in Fig. 12. Fig. 13 
indicates that the symmetry of the oscillation requires 
only a half-cycle to be studied. The boundary be- 
tween stability and instability of the solution is, then, 
replacing small deviations by differential deviations 


doy = —da, (46) 
Using Eq. (24) to trace differential deviations from 
(1) to (4), 


do, + dé,{1 — [1/(1 — 6)]} = doy + do {l — 
[1/(1 — o4)]} (47) 


Since sg < 5,, contours (1), (4), and (9) are all straight 
lines, with the slope 


do/ds = e# — 1 (48) 


determinable from Eq. (31). Since the trajectory be- 
tween contours (9) and (1) is also a straight line 


do, = day ; (49) 
do, = day (50) 


Substituting Eqs. (46), (48), (49), and (50) into Eq. 
(47), and solving for 6, 


G = 1 — e~4[(1 — o2)/(1 — 203)] (51) 
From Eq. (31) with a, b, k = 0 
= oi(e%4 > 1) (52) 


Combining Eq. (52) with Eq. (24) between contours 
(1) and (2), 


o,e°4 + In (1 > G}) = G2 + In (1 — G2) (53) 
Putting Eq. (51) into Eq. (53) 
ed — Sq = Ge 4+ In (1 -— 2a) + (1 on d2)/(1 cane 262) (54) 


The boundary defined by Eq. (54) has been plotted 
in Fig. 12. Above the boundary, the solutions given by 
the curves of constant o, represent modes of oscillation 
which the physical servo can attain and maintain. 

Below the stability boundary, the solutions represent 
boundaries themselves. If a servo is hunting with an 
amplitude slightly greater than that given by Fig. 12 
below the stability boundary, the amplitude will in- 
crease on successive swings, until the upper or “‘stable’”’ 
amplitude of oscillation indicated in Fig. 12 is attained. 
On the other hand, if the initial amplitude is slightly 
less than that given by Fig. 12 below the stability 
boundary, the amplitude will decrease on successive 
swings until the servo comes to rest. 

Designers of relay servos are familiar with the 
phenomenon of a relay servo that, when slightly dis- 
turbed from rest, will return to rest, but when given a 
vigorous impulse from rest, will hunt indefinitely in an 
oscillation of approximately constant amplitude. 








372 


For the servo elements from which Fig. 12 was de- 
duced, steady-state hunting can be prohibited for any 
magnitude of initial disturbance by making the contact 
spacing larger than the value given by the intersection 
of the stability boundary with the delay time value for 
that particular servo. 


STABILIZING DEVICES 


In order to transform the basic servo elements con- 
sidered in the preceding paragraphs into a high-per- 
formance relay servo, a stabilizing means more effective 
than simple contact spacing is required. The method 
selected must drastically reduce the amplitude of 
steady-state oscillation in the presence of delay time, 
particularly for the case of nonstationary input. 

Fry’ has described a number of stabilizing devices 
applicable to relay servos, These devices for damping 
the steady-state oscillation are conveniently divided 
into two classes, one class consisting of those devices 
whose description requires consideration of one or more 


additional degrees of freedom of the complete mecha- : 


nism, and the other class consisting of those devices 
whose operation is perfectly described by contours in 
the servo’s phase plane. An extension of the methods 
of the present study will handle the former class, in- 
cluding the Lanchester damper'* '® and the clutch- 
brake combination® for disengaging the load from the 
servo at zero error, but the extension will be reserved 
for a later paper. The latter group of stabilizers in- 
cludes the friction brake, flip-flop and front-lash de- 
vices, and anticipatory circuits. These will now be 
discussed in turn. 


Friction Brake 


For servos employed to match inputs that vary slowly 
or occasionally, a magnetically controlled friction brake 
may be employed to eliminate hunting about the zero 
position. The clutch coil indicated in Fig. la and here- 
tofore unmentioned controls such a brake. When both 
motor fields are open, the brake is applied to the load. 
When either field is energized, the brake is removed 
from the load. The actual servo from which Fig, la 
was drawn® also had a clutch controlled by the same 
circuit that disengaged the motor armature from the 
load at the instant the brake was applied to the load. 
If the duty of the servo is simply to move the load 
from one position to another temporarily constant, 
this sort of device is satisfactory in aiding the servo to 
bring its load to a “‘dead-beat”’ stop. 

When the input is varying rapidly enough so that the 
clutch circuit is continually being closed and opened, 
quantitative study of the motion requires that separate 
phase-plane plots of servo and load must be made across 
the open contact interval. Methods so far developed 
are not directly applicable to this case. 

If the inertia of the servo is small, compared with 
that of its load, so that no provision is made for un- 
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coupling the servo drive at zero error, usefulness of the 
friction brake can be studied by diagrams based on 
Fig. 3, in which the dry friction torque is augmented 
by the friction drag of the brake in the open contact 
interval. Lag contours can be computed for the study 
of the effect of delay time in torque control and brake 
operation, 

The friction brake applied to a high-performance re- 
lay servo has the disadvantage that it tries to lock the 
load to a fixed reference whenever it is energized. In 
order to reduce servo errors in following rapidly moving 
inputs, other stabilizing devices are preferable. 


Flip-Flop and Front-Lash Devices 


Referring to any of the phase-plane diagrams for 
hunting in the presence of delay time, it is evident that 
the steady-state oscillation follows clockwise rotation 
of the contours of torque reversal. The flip-flop and 
front-lash damping devices attempt to reduce the am- 
plitude of steady-state hunting by rotating the con- 
tours of torque reversal counterclockwise in the phase 
plane, 

The flip-flop circuit is similar in action to simple 
contact spacing. The difference is that the formerly 
inactive zone between the energized contacts is now ac- 
tivated in such a manner as to oppose the servo motion, 
regardless of the direction from which the servo ap- 
proaches this mid-segment. Fig. 14 shows how a flip- 
flop circuit may be fabricated for the servo of Fig. la 
by the addition of a relay. 

The flip-flop circuit does not prevent hunting but it 
does reduce its amplitude. Depending upon the rela- 
tive magnitudes of delay time and width of the mid- 
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Fics. 15a, 15b, and 15c. Modes of steady-state oscillation of a 
relay servo with “‘flip-flop’’ damper. 
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segment, the steady-state oscillation takes the form of 
one of three modes, The modes are shown in Figs. 
l5a,b, andc, For an extremely wide mid-segment, the 
mode of Fig, 15a appears, in which the servo simply 
oscillates about the boundary between the mid-segment 
and one of the outer segments. This will be called the 
“lower mode.” 
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For an extremely narrow mid-segment, the mode of 
Fig. 15c is performed, with both boundaries of the mid- 
segment crossed twice in each cycle. This will be called 
the “‘upper mode.” 

The intermediate mode, illustrated in Fig. 15b, has 
a number of submodes, depending upon how many 
cycles comprise the pattern. Successive cycles are not 
identical, as in the case of the other modes, but several 
cycles form a pattern that repeats indefinitely. 

A front-lash damping device is shown in Fig. 16. The 
point of torque transfer is advanced by moving the roll- 
ing contact relative to the contact plate with the aid of 
a dry-friction coupling. In Fig. 16 only the friction 
drag arm and the rolling contact touch the contact 
plate. The device is proportioned to advance the roll- 
ing contact at k; times the servo rate, up to the limits 
of the stops. 

The front-lash stabilizer permits three modes of os 
cillation, These are illustrated in Figs, 17a, b, and c. 
The upper mode, shown in Fig. 17c, is identical to the 
upper mode for the flip-flop stabilizer, when the stop 
spacing o, equals the mid-segment width o,,. 

The lowest mode, shown in Fig. 17a, is performed 
when the amplitude of oscillation is so small that the 
front-lash arm touches neither stop. The case then re- 
duces to that of an elementary relay servo with delay 
time and zero cqntact spacing, and the amplitude may 
be computed accordingly. 

The intermediate mode is obtained when contact 
transfer takes place before the front-lash arm has 
crossed the full-stop width, but the remaining travel 
takes place before the servo velocity is reduced to zero, 
This mode is similar to the upper mode except that, in 
the case of the latter, the front-lash arm has crossed 
the full-stop width before contact transfer takes 
place. 


Mathematical Solutions 


Figs. 15 and 17 are drawn for the case of stationary 
input and constant delay time in torque control, The 
methods used to determine the amplitude of oscillation 
of the relay servo without stabilizer can be applied di- 
rectly to obtain the amplitude with stabilizer. The 
upper mode of oscillation is the same for both flip- 
flop and front-lash dampers, and will be considered 
first. 

Assuming that a steady-state oscillation exists in 
this mode, and solving for its amplitude, this is found 
to be given by 


erd[e— me] — Go)] — 2e4 + (1 — G2) =O (55) 


Here o,, and o, are used interchangeably, depending 
upon the stabilizer under consideration. Stability of 
the solutions indicated by Eq. (55) is determined by 
tracing a differentially displaced trajectory through 
half a cycle, and obtaining the boundary equation as 


(56) 


e~*4 = (1 — 2¢2)/(1 — 2)? 








374 JOURNAL OF THE 


CURVES OF CONSTANT CO, 
Pi CURVES OF CONSTANT CO, 










b 





oo 
$4 ——> 


Fic. 18. Amplitude of steady-state oscillation of a relay servo 
with inertia, viscous damping, constant delay time, ‘‘flip-flop” 
damper, and stationary input. 


Solutions of Eq. (55) for the relay servo with flip- 
flop stabilizer have been plotted in Fig. 18. The sta- 
bility bound is also shown. It has the same significance 
in this case as it had for Fig. 12, with the exception that 
in Fig. 18 it enables the reader to determine when the 
oscillation deteriorates to a lower mode. 

No solutions will be carried out for the lower modes 
for the flip-flop stabilizer. Since the torque advance 
occurs less than twice in each cycle in the lower modes, it 
will be assumed that the stabilizer is correspondingly 
less efficient. 

The boundaries of the three modes of oscillation of 
the relay servo with front-lash stabilizer can be written 
down by inspection of Fig. 17.. They are 


o, < —kiog Upper Mode (57) 
—kyo, < o, < —2k,o, Intermediate Mode (58) 
o, > —2k\c, Lower Mode (59) 


From Fig. 15c it can be observed that the solutions of 
Fig. 18 for the upper mode are also applicable to the 
intermediate mode of oscillation of the front-lash sta- 
bilizer, if the stop width a, is replaced by an effective 
stop width (oz)eg., given by 

(er)ea. = —2hioq — O71 (60) 

For the particular case of k; = —1.00, Fig. 18 has 
been modified to indicate the amplitude of oscillation in 
all three modes of a relay servo with front-lash stabili- 
zation. Note that Fig. 19, displaying the results, in- 
dicates that for a given value of delay time sy, there 
is an optimum value of stop spacing oz, occurring 
at the lower boundary of the upper mode of oscilla- 
tion. 

Although the case has not been studied in detail, 
when the amplitude indicated by Fig. 19 lies below the 
stability boundary of Fig. 18, it is likely that a mode of 
oscillation appears which combines the characteristics 
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Fic. 19. Amplitude of steady-state oscillation of a relay servo 
with inertia, viscous damping, constant delay time, ‘‘front-lash” 
damper, and stationary input. 


of Fig. 15b and Fig. 17b. One may reasonably expect, 
from considerations similar to those applied to the flip- 
flop stabilizer, that the optimum values of &; and oz, for 
a particular delay time, are those which place the am- 
plitude of oscillation in Fig. 19 exactly on the stability 
bound. This can be ascertained graphically, but has 
not been carried out at the time of writing. 


LEAD STABILIZATION 


Under what has been assumed to be the most favor- 
able condition of their operation—hunting at the lower 
bound of the upper mode of oscillation—the front-lash 
and flip-flop dampers operate by advancing the contour 
of contact interchange horizontally along the phase 
plane’s o axis. For a constant delay time, however, 
the servo travels a‘ greater distance between contact 
transfer and torque transfer at high velocities than at 
low. The contact transfer contour should therefore be 
advanced farther ahead of the ¢ axis for large velocities 
than for small. 

The effect can be produced by a generator driven by 
the servo.’ The contacts of Fig. 2 are replaced by a 
potentiometer, or similar device, displaying a voltage 
proportional to the error. The relay of Fig. la is 
controlled by the difference between voltage propor- 
tional to error, and voltage proportional to servo veloc- 
ity, so that relay operation takes place when 


Tjbn — (0m — 9.) = 0 (61) 


where 7, is a positive constant determined by the cit- 
cuit elements.and subject to control by the designer. 
Eq. (61) defines a straight line contour of relay opera- 
tion in the phase plane, leading the ¢ axis by an angle 
determined by 7,. A purely mechanical combination 
of spring plus viscous damper can also be employed to 
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utilize Eq. (61) without the generator and associated 
additional electrical equipment.’ 

Since the contour defined by Eq. (61) is a straight 
jine, Eq. (31) still defines the straight line of actual 
torque reversal when there is a constant delay time. 
Graphical construction of trajectory webs for represen- 
tative relay servos with lead stabilization and delay 
time quickly reveals that the stabilization defined by 
Eq. (61) causes the servo to follow a constant velocity 
input with a steady-state error whose mean is roughly 


5:5 = wT, (62) 


One way of eliminating this steady-state error is to 
cause the input itself to drive a generator, independent 
of the servo-driven generator, and to control the relay 
by the difference voltage between the two units. When 
the input has not sufficient torque to drive a generator, 
the same result can be obtained by passing the voltage 
proportional to error through a differentiating circuit, 
then a vacuum tube amplifier if necessary,® so that the 
relay of Fig. la is operated when 


7,6 +6 =0 (63) 


This system has no steady-state lag in following a 
constant velocity input. 


MATHEMATICAL SOLUTION 


For the particular case of stationary input, the 
mathematical determination of amplitude of oscilla- 
tion of a relay servo with lead stabilization is easily 
carried out. For this case Eqs. (61) and (63) are iden- 
tical and results apply to both types of stabilizers. 
Writing Eq. (63) in dimensionless form 


so +o =0 (64) 


where s, = T,/r. 
The amplitude of steady-state oscillation is deter- 
mined by the methods earlier described as 


e*4[(1 = oe? 24-0) | d¢% +(1—o.)=0 (65) 


Fig. 20 shows solutions of Eq. (65). The lead sta- 
bilizer does not eliminate steady-state hunting but will 
reduce its amplitude to as small a value as desired, if 
s, is chosen large enough. The fact that hunting is 
not eliminated completely is a consequence of the as- 
sumption of constant delay time in torque control. As 
long as the quarter period of oscillation is large com- 
pared with the delay time, this is an adequate assump- 
tion, When the magnitudes are reversed, however, 
the exact character of the delay becomes of importance 
in obtaining a true representation of the physical sys- 
tem. While the operation of a relay probably involves 
a fairly constant delay time, torque generation in a 
d.c, series motor may be better represented as 


T(t) = Tax.(1 — €~”™") (66) 


where ¢ is measured from the instant of contact closure, 
and 7,, depends upon the inductance-resistance values 
of the motor windings. 

Eq. (66) represents a type of lag which is less bother- 
some than constant delay time, since some torque ap- 
pears almost immediately after contact closure. One 
might consider the full torque as divided into a number 
of parts, with some of which there is practically no delay 
time associated. This is a more favorable condition 
than that of constant delay time, where the delay time 
is associated with the total torque. When the physical 
servo is subject to lag that tends to follow the law of 
Eq. (66), it may be expected that computations per- 
formed by the methods of this paper for constant delay 
time will be highly conservative, for high-frequency 
oscillations. The result may be complete suppression 
of the oscillations which, in Fig. 20, are merely reduced 
to negligible amplitudes. 

If, however, the lag is truly a constant delay, steady- 
state oscillation can still be suppressed by employing 
a more complex network and amplifier to replace Eq. 
(63) by 


f(o) +o =0 (67) 


It may be deduced from sketches in the phase plane, 
that a suitable f(c) is 


f(a) = ket + ks sgn (c) (68) 


where ko, kz are constants. 


COMPARISON OF DAMPERS 


It is of considerable interest to compare the effective- 
ness of some of the stabilizing devices in reducing the 
amplitude of steady-state oscillation. The case of 
stationary input has been selected, because of the avail- 
able mathematical solutions, but only those stabilizing 
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Fic. 20. Amplitude of steady-state oscillation of a relay servo 
with inertia, viscous damping, constant delay time, ‘‘lead’’ sta- 
bilization, and stationary input. 
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Fic. 21. Comparison of the effectiveness of several types of 


dampers in suppressing steady-state oscillation of a relay servo. 


devices have been chosen which can be applied to servos 
with rapidly moving inputs. Fig. 21 shows the results. 

Ordinates are the half-amplitude of oscillation, or an 
equivalent measure. For example, Fig. 12 reveals that 
in the case of simple contact spacing, the total amplitude 
of oscillation decreases uniformly as contact spacing is 
increased, until the spacing is large enough to stop the 
oscillation. The uncertainty in servo position is then 
the value of contact spacing, and half this value has 
been plotted in Fig. 21. 

For the flip-flop stabilizer, it is observed from Fig. 
18 that the total amplitude of oscillation decreases uni- 
formly in the upper mode as the width of the mid-seg- 
ment is increased, until the stability bound is reached. 
Below this bound a mode of oscillation appears in which 
the stabilizing action takes place less frequently in each 
cycle and, since it is known that when the lowest mode 
is reached the uncertainty in servo position is greater 
than that of the ‘“‘no-stabilizer’”’ curve in Fig. 21, the 
amplitudes corresponding to the stability boundary in 
Fig. 18 have been accepted as indicative of the perform- 
ance of the optimum flip-flop stabilizer. Similar 
reasoning leads to the selection of the same curve for 
the frontlash stabilizer. The half-amplitude plotted 
in Fig. 21 is 


1/, = Og + 1/ 90m (69) 


The lead stabilizers are assumed to have s, a constant 
ratio to sz, on the observation that large anticipation is 
necessary only when there is large delay time. 

The diagram shows the sort of results that might 
have been expected, The flip-flop stabilizer, which ap- 
plies an opposing torque across its mid-segment, is 
about twice as effective as contact spacing, which 
merely removes the driving torque over the interval. 
The flip-flop stabilizer is able to reduce the amplitude 


of steady-state hunting to 40 per cent of its magnitude 
with no stabilizer. 

The lead stabilizer, however, permits the amplitude 
of oscillation to be reduced to any value desired. It is 
the most logical damping device to be applied to a high- 
performance relay servo. 


CONCLUSION 


The methods developed in this paper have the con- 
siderable advantage that they subject to quantitative 
analysis many factors of servomechanism design which 
could not be studied heretofore by simple processes, 
The flexibility of the phase-plane approach, particularly 
in its graphical form, is so great that it has been pos- 
sible to explore only a few of its applications in the pres- 
ent paper. 
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An Analogy Among Torsional Rigidity, 
Rotating Fluid Inertia, and Self-Induc- 
tance for an Infinite Cylinder 


JOHN W. MILES* 
University of California at Los Angeles 


SUMMARY 


It is shown that the cross-sectional, polar moment of inertia 
(I) of any cross section containing a perfect fluid is given by 
subtracting froth the polar moment of inertia of the cross section 
(I4) (as a rigid body) the torsional rigidity of the cross section 
(J), the latter being calculated as though the cross section were 
perfectly elastic with a shear modulus of unity. Results are 
given for an ellipse, an equilateral triangle, an isosceles right 
triangle, and a rectangle. References are included from which 
results for various other shapes may be obtained. 

It is also shown that the quantity J is given by the self-induc- 
tance of a conducting cylinder of high permeability at low fre- 
quencies. No applications are made of this last analogy. How- 
ever, it is felt that it might often be more convenient than the 
membrane analogy for determining torsional rigidities. 


INTRODUCTION 


I CARRYING oOuT flutter and stability calculations, 
it is often desirable to know the inertia of a wing 
containing fuel. The added inertia of the fuel in trans- 
lation for a full tank may be taken as if the fuel were 
solid, but, inasmuch as normal fuels are characterized 
by small viscosity, the rotational inertia is smaller 
than it would be if the fuel were solid. There follows 
a theoretical calculation of the rotational inertia for 
a perfect fluid. 


THE STREAM FUNCTION 


The problem of obtaining the stream function for 
a rotating cylinder containing a liquid has been con- 
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sidered by Milne-Thomson.' If (u,v) are the x and y 
velocity components of a fluid in the (x,y) plane, the 
stream function y is defined by 


u = —dy/dy (la) 
dy / dx (1b) 


v 


As is well known, the stream function defined in Eqs. 
(1a) and (1b) satisfies Laplace’s equation 


Vy =0 (2) 
For a uniform rotation of angular velocity w about 


the origin, it is easily seen’ that the stream function 
is 


Vs = (w/2)(x? + y?) + A (3) 


where A is an arbitrary constant.t For the problem 
at hand, Eq. (3) is the boundary condition imposed on 
y by the cylinder. It is then expedient to choose y 
in the fluid so that 


¥(x,y) = V.(%,y) + 2wx(x,y) (4) 
where 
x(x,y) = 0 (5) 
on boundary. 
Substituting Eqs. (3) and (4) in Eq. (2), the partial 
differential equation for x is seen to be 


V*x(x%,y) = —1 (6) 

t For a cross section bounded by a singly connected curve, 

it is convenient to choose A = 0, although this may not be so 
for a multiply connected curve 


From Eq. (1) it may be seen that the square of the magnitude of the velocity (u* + v?) is given by the gradient 
squared of the stream function; hence, separating y into y, and x, the kinetic energy (7) per unit length of 


cylinder, having the cross section S, is given by 


T = YoalS S3(Vvu)*dS + 40S S5V vo'V xdS + 40° SL S(V x)*d 5] (7) 


The first integral in Eq. (7) is the kinetic energy calculated as if the fluid were a solid moving with the boundary 


and gives the corresponding moment of inertia (Z,). 


The second and third integrals may be modified by apply- 


ing Green’s first theorem” ‘ to an infinitesimally thin lamina of the cylinder, viz: 


S SN ovV dS = —S SobiV *bxdS — $ o1(H-V $2)dl (8) 
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where 7 is a unit vector normal to the boundary of the cylinder. The second integral in Eq. (7) may be evaluated 
by letting ¢; = x and ¢: = y,, while the third may be handled by putting ¢; = ¢. = x. Observing, from Eq. 
(5), that x vanishes on the boundary and, from Eq. (3), that V*y,, = 2w, while, from Eq. (6), V?x = —1, Eq. 
(7) yields 


I = 2T/po? = 1g — 4S foxdS = 1, — J (9) 


where J is the cross-sectional, polar moment of inertia about the origin and J, is defined above. 

The problem of evaluating the effective inertia of the (perfect) fluid in a rotating cylinder has been reduced to 
that of obtaining a solution to Eq. (6), satisfying the boundary condition, Eq. (5), and integrating it over the 
cross section of the cylinder. Eq. (6) is Poisson’s equation for a uniform source distribution and arises in such 
familiar problems as the torsion of a cylindrical bar,? the deflection of a homogeneous membrane subjected to 
uniform loading,? and the flow of a uniform current through a cylindrical wire.* Hence, for those cross sections 
that do not lend themselves to simple analysis, the problem at hand may be solved by analogy with problems in 
other fields which have previously been solved, either analytically or experimentally. Two analogous problems 





will be treated below. 





TORSION OF CYLINDRICAL BARS 


The problem of torsion in a cylindrical bar, as first 
solved by Saint-Venant, has been treated by Timo- 
shenko.*? Briefly, if ¢ is defined as a ‘‘stream function”’ 
for the shear stress equation (T132),* it obeys the 
equation (T133) 


V%=F (10) 


where F is an undetermined constant. Moreover, 
from (T134), @ may be taken to vanish on the boundary 
of a cross section bounded by a singly connected curve. 
The twisting moment is given by equation (T135) as 


M = 2S S5odS (11) 
while the angle of twist, from equation (T139), is 
6 = —F/2G (12) 


where G is the shear modulus. Setting ¢ = — Fx and 
substituting in Eqs. (10) through (12), the torsional 
rigidity of the cross section is given by 


C= M/@=4GS fyxdS = GJ (13) 


where x is determined by Eqs (5) and (6). Hence, 
J in Eq. (9) is given by the torsional rigidity of the 
cross section for a shear modulus of unity. Certain 
applications of this analogy have been discussed by 
Timoshenko? who attributes the analogy to Greenhill.* 


INTERNAL INDUCTANCE OF A CYLINDER 


Consider a cylindrical conductor of cross section S 
and (absolute) permeability u, carrying an axial current 
of sufficiently low frequency to permit the assumption 
of uniform current density over S. The magnetic 
field components (Bx, By) may then be computed from 
a stream function (A) which obeys the equation (using 
MKS units) 


VA = —pi (14) 


* The equations preceded by a “T”’ are from reference 2. 


‘permeability, being dimensionless. 


where 7 is the current density.‘ (A is the z component 
of the vector potential 4, whose curl is the magnetic 
field B.) The energy of the magnetic field per unit 
length within the conductor is given by* 


W = (1/2u) Sf fGB°dS = pf SiAdS = "/oL,,(iS)* 
(15) 


where L, is defined as the internal self-inductance of 
the conductor in henries, and (2S) is the current in 
amperes. Substituting x = (A/i) in Eqs. (14) and 
(15) and setting » = 4% X 10-*(u/uo), where yo: is the 
permeability of free space, L; is given by 


+ Ly = (p/S)S SgxdS = [10% (u/u0)S-*]J (16) 


Hence, the integral J is given by multiplying the 
internal self-inductance (henries) by 10°r~'(u/yo) ~1S?, 
S being in square meters and (u/yo), the relative 
Comparing Egs. 
(6) and (14) it is seen that x satisfies Eq. (6). In order 
to yield the corresponding boundary condition, Eq. 
(5), it is necessary that A vanish (or be constant) along 
the boundary of S, which implies that the magnetic 
field is tangential at the boundary. The latter condi- 
tion will be approximately satisfied for a conductor of 
high permeability; moreover, for high permeability 
the internal self-inductance is practically the total 
self-inductance, since the field outside the conductor 
is negligible. 

The author doubts that this analogy will prove 
useful in determining J, but it is often desirable to 
determine L, for certain odd cross sections at low 
frequencies for transient analyses, in which case J, 
might well be determined by a torsional measurement 
of J. It is hoped to carry out experimental measure- 
ments at some future date. 

The so-called ‘‘penetration depth” of any conductor 
(MKS units) is given by 6 = (ufo)—/*/em (reference 4), 
where f is the frequency, o is the conductivity, and pis 
the absolute permeability. For iron of relative per- 
meability 10,000, this gives an approximate depth 
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Fic. 1 (top left\ Ellipse. 
triangle. Fic. 3 (bottom left). 
(bottom right). Rectangle. 
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Fic. 2 (top right). 
Right isosceles triangle. 


Equilateral 
Fic. 4 


(meters) of 5 = 1.58 X 10-3f-"/, and it is clear that 
the foregoing analogy will be useful only for direct cur- 
rents. 


EXAMPLES 


In order to solve Eqs. (5) and (6), it suggests itself 
to choose 


x(x,y) = k(x,y) f(x,y) (17) 


where f(x,y) = 0 is the equation of the contour of S, 
since such a function automatically satisfies Eq. (5). 
Thus for the ellipse shown in Fig. 1: 


S(x,y) = [(%/a)? + (y/b)? — 1] = 0 (18) 


k, as determined by substitution in Eq. (6), is given by 
the constant 


k = —1/20°*(a? + 


Substituting x, as given by Eqs. (17) through (19), 
in Eq. (9) yields 


J = ra*h*(a? + 6?)—! (20) 


For a circle of radius a(= 0), Eq. (20) reduces to J = 
ra‘/2, which is the polar moment of inertia of the cross 
section. 

The simple case of a circle presents a convenient 
check on the foregoing analogies. For a rotating 
circular tank, Eq. (9) gives an effective fluid inertia of 
zero, and it is physically obvious that this result is 
correct for zero viscosity. It is well known? that the 
torsional rigidity of a circular shaft is simply its moment 
of inertia (multiplied by G). The internal self-induct- 
ance of a circular wire of permeability u is given by 
0.5 X 10-(u/uo) henries;? substituting this result in 


b*)~" (19) 


Eq. (16), together with S = xa’, yields J = ra‘/2, in 
agreement with Eq. (20). 

Other methods of obtaining solutions to Eqs. (5) 
and (6) have been discussed extensively in the litera- 
ture’ ** and will not be dealt with here. Timo- 
shenko’ derives results for many useful cases. For 


the equilateral triangle of altitude a shown in Fig. 2 


J = a*/15v/3 = 0.61, (21) 


where J, is the polar moment of inertia about the 
centroid. For the right isoSceles triangle of small side 
a shown in Fig. 3 


= a*/38.3 = 0.4701, (22) 


where J, is again referred to the centroid. For the 
rectangle of large side 2b and small side 2a shown in 
Fig. 4 





= k,(2a)*(2b) (23) 
mb 
1 192 = — 
b = =) 1- — (2) 5: senadnantin (5) (24) 
3 mw \b/ n=0 
(2n + 1)5 
For a narrow rectangle (b >> a), Eq. (23) reduces to 
J = '/;(2a)3(2b)[1 — 0.630(a/b)) (25) 
while for a square (a = )) 
J = 0.406(2a)* = 0.844I, (26) 


A plot of 3k, as given by Timoshenko, is given in Fig. 5. 

Additional results are given by Timoshenko for a 
sector of a circle, several cross sections defined by 
f(x,y) = 0 where f(x,y) is a polynomial, and for thin- 
walled tubes. Application of the methods of Ritz 
and Trefftz to the determination of x(i.e., ¢) is also 
discussed. 

Returning to the original problem, the effective 
moment of inertia of the fluid in a rotating elliptic 
tank is given by substituting 











I, = (x/4)ab(a? + 5?) (27) 
and J from Eq. (20) in Eq. (9) to obtain 
(I/Ia4) = [(@? — 6*)/(@* + 5%)}* = * (28) 
3K, 
10 J 
os] 
061 
O* + t t + " 
o aA, 


02 o-4 06 os 
Fic. 5. XK, for a rectangle. 
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(It may be observed that problem 12 on page 251 of 
reference 1 gives the result, Eq. (28), with the quantity 
in brackets cubed, probably a typographical error.*) 

For most practical purposes the effective inertia of 
a rotating cylinder of fluid can probably be determined 
with sufficient (i.e., warrantable) accuracy by replacing 
the cross section by an equivalent ellipse or rectangle 
and applying the foregoing results. An equivalent 
ellipse should be particularly convenient because of the 
simplicity of Eq. (28). Moreover, ‘Saint-Venant 
found that the torsional rigidity can be calculated 
approximately. ...by replacing the given shaft by the 
shaft of an elliptic cross section having the same cross- 
sectional area and the same polar moment of inertia 
as the given shaft bar.’’? 

The application of the foregoing to fuel tanks is, of 
course, limited by the assumption of a perfect fluid. 
However, normal fuels possess small viscosity, and the 
results should be satisfactory, provided that the tank 
is sufficiently full to prevent a great amount of splash- 
ing. (This latter requirement is, of course, better 


* It has been called to the author’s attention that this result 
is stated correctly on p. 87 of Lamb’s Hydrodynamics. 


satisfied for small amplitudes of motion and high fre. 
quencies, where the effects of gravity are negligible) 
If the tank is only slightly empty, the resultant fre 
surface implies that the effect of viscosity at the wall 
will be less important and the conditions of a perfec 
fluid are more closely approached. The justificatioy 
of applying the results obtained for a uniform rotation 
to nonuniform motion, as in the case of flutter, follows 
directly from the assumption of an incompressible 
fluid. 

It should be emphasized that Eq. (9) applies fo; 
any axis of rotation, and, although J, depends on the 
location of this axis, J depends only on the shape of 
the cross section. 
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Letter to the Editor 


(Continued from page 363) 


It can be stated generally that the buckling strain of a column 
is independent of its material but only in the elastic range. The 
physical reason for the multiplier (Ereg./Esec.) in the inelastic 
range should be clear from a strain energy consideration of the 
buckling phenomenon. As is known, the buckling load times the 
relative displacement of the end points of the column due to flex- 
ure at buckling is equal to the strain energy stored in, the column 
while it undergoes deflections corresponding to the buckled shape. 
When buckling occurs above the elastic limit, the flexural rigidity 
of the column is E,-q.J rather than EJ, where Erea. is the “aver- 
age’”’ between Young’s modulus £ and the tangent modulus E,, 
calculated from the Engesser-K4rman formula Ereg. = 4EE,/- 
(VE + VE,)?. (The formula is rigorously correct only for rec- 
tangular sections, but it holds approximately for other sections 
also.) 

The end load acting upon the column at buckling is o,.A, 
where A is the cross-sectional area of the column. Solution of 
the strain energy equation, expressed earlier in words rather than 
symbols, yields Eq. (6). In the gradual process of loading the 
stress and strain in the straight column increased in accordance 
with the stress-strain curve of the material. Hence the critical 
strain ¢,.,. can be obtained from the critical stress given in Eq. (6) 
by dividing the latter by the secant modulus Exec... The division 
yields Eq. (5). 

Although the secant modulus method for determining plate 
instability thus lacks theoretical foundation, it is perfectly ac- 
ceptable in the case of materials for which its validity is estab- 
lished by sufficient experimental evidence. It should be recalled 
that the validity of Eq. (6)—and consequently of Eq. (5)—was 


first demonstrated in the experiments carried out by von K4rmén 
and published in 1910. Although his results have been corrob- 
orated by several other investigators, it is common practice to- 
day in the airplane industry to use the tangent modulus rather 
than the reduced modulus in the Euler formula. This practice is 
conservative and hence allows for inaccuracies of manufacture, 
If an empirical formula is preferred in column analysis where a 
perfectly reliable theoretical formula has been available for almost 
40 years, there is every reason to accept empirical results ex- 
pressed, for instance, by the secant modulus method in plate 
buckling calculations where an equally well-founded theory is 
lacking beyond the elastic limit. 

As far as the writer’s suggestions regarding the calculation of 
the buckling stress of plates are concerned (‘A Note on Inelastic 
Buckling,” published in April, 1944, in the JouRNAL OF THE 
AERONAUTICAL SCIENCES), they were only intended for use when 
a less conservative solution is not conveniently available. The 
writer is willing to drop his suggestion in favor of any simple em- 
pirical relationship that is sufficiently corroborated by test. It 
is Mr. Gerard’s contribution to have attacked by an elegant ex- 
perimental method the somewhat neglected though important 
problem of plate buckling beyond the elastic limit, to have sug- 
gested the simple solution of using the secant modulus in the for- 
mulas, and to have justified his suggestion at least in the case of 
24ST aluminum alloy by the good agreement found between his 
method and the results of his experiments. 

. N. J. Horr 
Associate Professor of Aeronautical Engineering 
Polytechnic Institute of Brooklya 
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Dynamic Loads in Airplane Structures 
During Landing | 


SHOU-NGO TU* 
Curtiss-Wright Corporation 


SUMMARY 


The additional dynamic loads introduced by vibrations excited 
by the impact load in landing may create a critical design condi- 
tion for the structures of an airplane. In this paper, a method is 
developed to estimate the magnitude of these loads in the wing 
for the two-wheel landing condition. The wing is found to 
vibrate in such a manner that the points of attachment of the 
landing gear are always two of the node points. This fact gives a 
guide to determine the coupled modes and frequencies. For- 
mulas are derived then toshowtherelativecontributions of various 
modes to the resultant motion. By calculating the inertia loads 
and moments from the accelerations thus determined, the beam 
shear and bending moment can be easily evaluated. 


NOTATIONS 


distance of the point of attachment of the landing 
gear from the plane of symmetry of the wing, in 
ft. 

an element in the frequency determinant 


a; = 

Ay = the cofactor of a;; in the frequency determinant 

B = bending rigidity of the wing, in Ibs.ft.? 

3 = torsional rigidity of the wing, in Ibs.ft.? 

I = moment of inertia about the line joining the points 
of attachment of the landing gear to the wing, 
in slug-ft? 

K,(é) = a function denotes the relative magnitude of the 
nth coupled mode, Eq. (16) 

L = wing semispan, in ft. 

m = mass, in slugs 

8 = the order of vibration mode 

P = inertia force, in Ibs. 

P(a) = reaction at the landing gear due to the inertia 
forces, in Ibs. 

| = integer between 0 and r, denoting the location 
of masses in the wing 

?.4q = distance from plane of symmetry, in.ft. (Used 
in the Appendix only) 

r = number of degrees of freedom used in the vibra- 
tion analysis 

R(p) = slope at the support of the beam due to a unit load 
at station p : 

s = the operator 0/0 

S = static moment of inertia, in slug-ft. 

t time, in sec. 

§ inertia twisting torque, in lbs.ft. 

U,V quantities defined by Eqs. (12) 


V(0) 
Vii) 


airplane landing speed, in ft. per sec. 

velocity of the point of attachment of the land- 
ing gear 

velocity of the wing 

the Laplace transformation of V(t) 

distance from plane of symmetry of the wing 


V(x, t) 
V,(s) 


x 
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vertical displacement of the beam, in ft. 


y = 
Y(x,w,) = vertical displacement in the nth coupled mode 
Yuo(a, w,) = the change of Y with respect to w, Eq. (14) 
Yi(x,s) = the Laplace transformation of y 

T = the parameter for time 

r = torsional influence coefficient 

bb = bending influence coefficient 

Wr = the vibration frequency of the mth coupled mode 
0 = angular displacement of the beam, in rad. 
O(x,w,) = angular displacement in the mth coupled mode 
@,(x, s) = the Laplace transformation of @ 


INTRODUCTION 


AN AIRPLANE LANDS, the vertical component of 
its velocity is reduced to zero in a short interval of 
time. This deceleration introduces inertia loads 1n the 
structures. In the past it was usually assumed that the 
inertia load factor was the same for all parts of the air- 
plane. This assumption is approximately true if the 
size of the airplane is small and the structures are com- 
paratively rigid. For larger airplanes, especially those 
with engines installed in the wing, vibrations may some- 
times be excited by the impact force transmitted 
through the landing gear. These vibrations will intro- 
duce additional inertia loads which, in some instances, 
may create a critical design condition for the structure. 
In this paper a method is presented to estimate the 
magnitude of these loads in the wing. In a similar 
manner, the inertia loads in the fuselage and tail can be 
calculated, but the details are not discussed here. 


A 


THE LANDING PROBLEM 


The airplane under consideration is a conventional 
landplane, with main landing gears directly attached to 
the wing structure. To simplify the analysis, the pres- 
ence of the nose or tail wheel is ignored. This simplifi- 
cation will not involve too many restrictions, since very 
often the auxiliary landing gear does not touch the 
ground for a short interval after the airplane lands. 
Even in a three-point landing, because of the longer 
distance from the c.g. of the airplane, the force trans- 
mitted through the auxiliary gear is comparatively 
small. Its neglection is therefore justified at least for 
the calculation of inertia loads in the wing. 

Because of the nonlinear characteristics of the landing 
gear, it is difficult to treat the whole airplane as a unit 
in the landing:problem. This difficulty can be avoided, 
however, by separating the dynamic characteristics of 
the landing gear from that of the elastic structure. 
The landing gear can be regarded as a damped spring; 
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the airplane structure, as an elastic beam with dis- 
tributed masses and moments of inertia. The beam 
is supported by the springs at two points, P and P’ in 
Fig. 1, symmetrically located from the midpoint of the 
beam. Before the wheels touch the ground, both the 
springs and the beam are descending at a constant 
velocity V(0). Upon touching the ground, the velocity 
of the springs is reduced to zero in a certain interval; 
but at the same time, the beam moves at a velocity 
V(x, t). The nature of V(x, ¢) depends upon (1) the 
characteristics of the landing gear, (2) the elastic prop- 
erties of the airplane structure, and (3) the initial 
velocity V(0) and the manner in which the airplane 
lands. Once V(x, ¢) or its equivalent is determined, 
the distribution of inertia loads will follow immediately. 
































ae 


Wing as a beam supported at points of attachment of 
landing gear. 


Fic. 1. 


Among the three items that determine V(x, #), the 
initial velocity V(0) is given. The elastic properties of 
the structure can be calculated in the usual manner. 
The dynamic characteristics of the landing gear can be 
represented by functions V(a, 4) and V(—a, #), the 
velocity of the points of attachment of the landing gear 
to the wing. V(a, t) does not necessarily equal V(—a, 
t); their relative magnitudes depend upon the manner 
in which the airplane lands. Only in symmetrical 
landing will their values be the same at all times. In 
spite of this difference, these two functions have the 
same fundamental characteristics. Mathematically, 
they can be expressed by the same function, with a 
proper phase difference and amplitude ratio. The form 
of that function can be determined by a drop test of the 
landing gear in question. For the purpose of this paper, 
both V(a, ¢) and V(—a, ¢) are assumed to be known. 


1946 


Since V(a, t) and V(—a, #) are, in general, complicated 
functions of time, it is not convenient to use them 
directly in the analysis. As Hooke’s Law is expected to 
hold in the elastic structures, the principle of super. 
position can be applied here. Therefore V(a, t) and 
V(—a, t) can first be decomposed into symmetrical and 
antisymmetrical components. The effects of each 
component are examined separately and final answer js 
obtained by superposition. 


THE CALCULATION OF V(x, ¢) 


As proposed in the preceding section, the airplane 
structure is replaced by an elastic beam, with its axis 
passing through the points of attachment of the landing 
gear to the wing. When the functions V(a, ¢) and 
V(—a, t) are given, the landing problem is reduced to 
finding the manner in which the beam vibrates when the 
motion of two of its points, P and P’ in Fig. 1, are 
prescribed. Since the beam has a plane of symmetry, 
only half of it needs to be considered under either sym- 
metric or antisymmetric loading. Taking the origin 
at the plane of symmetry and the x-axis along the 
beam axis, the dynamic equations for the beam in 
vibration are: 


0? °2) _ OY 079 
Ox? ( Ox? st ‘ies ri 


or? ot? 
(1) 
ra) 00 070 o? 
(Bre 9% 


at as” wre" 

where B and C are the bending and torsion rigidities of 
the beam (wing), m is the mass distribution, S and J are 
the first and second moment of inertia about the beam 
axis, y is the vertical displacement of the beam axis, and 
6 is the angular rotation about the beam axis. In the 
equations the damping is neglected. This is justified 
because the maximunt inertia load will occur in a 
fraction of a second and the damping effects will not be 
significant. The boundary conditions of the problem 
are: atx = L 


0*y/Ox? = 08y/dx* = 0 
0, under symmetric bending, 
dy/Ox = 08y/Ox® = 00/0x = 0 
at x = 0, under antisymmetric bending, 


y = 0*y/dx? = 0 = 0 


atx 


atx =a, 
Oy o*y 
VY Bx? Dax?” 
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dy/dt = V(t) 


at t = 0, 


y = 6 = 00/d0t = 0 
Oy/Ot = V(0) 


and for the beam as a whole, 


o°6 dy 
fe or + S32) ~=9 


where V(t) represents either the symmetrical or the 
antisymmetrical component of V(a, /). This problem 
can be solved more easily with the help of operational 
calculus. Let Yi(x, s) and 0,(x, s) be the Laplace 
transformation of y(x, ¢#) and @(x, #), respectively;' 
Egs. (1) are transformed into 


d? d?y, 
—{B 
dx? (8% 


d (Cde 
—_ (S ') + Is?0, + Ss?¥, = SV(0) 
dx \ dx 





2) + ms? Y, + Ss*O,; = mV(0) 
(2) 





All boundary conditions not involving derivatives of 
time are transformed into corresponding conditions by 
replacing partial differentiation with ordinary differ- 
entiation, y with Y,, and @ with 0,. The boundary 
condition that involves time is transformed into: at 
r=a 


sY,(a, s) = Vi(s) (3) 


where V;(s) is the Laplace transformation of V(t). 
Eqs. (2) can always be solved. Let the solutions that 
satisfy all boundary conditions but that represented by 
Eq. (3) be 

Vi(x, s) = KY(x, s) + [V(O)/s?] 


@1(x, s) = KO(x, s) 


(4) 


Where K is a constant still undetermined, Y(x, s) and 
0(x, s) are the solutions of the homogeneous simultane- 
ous differential equations corresponding to Eqs. (2). 
To satisfy the condition represented by Eq. (3), K is 
found to be 


C = {[Vi(s)/s] — [V(0)/s*]}[1/¥G, s)] 
Hence, from Eqs. (4), Y; and 0, become 





é V . s 7 
Yi(x,s) = Ee ne | . ae . he 
Vi(s) V(0)7 O(x, s) (5) 
iS <3 
saneeiit ow a, | Y(a, s) 


The functions y(x, ¢) and 6(x, ¢) can then be obtained 
from Eqs. (5) by summing up the residues. It is 
noticed that both Y; and ©, have simple poles at values 
of s which satisfies the equation Y(a, s) = 0, But at 
s = 0, Vi(x, s) has a pole of the second order, while 
(x, s) has only a simple pole. This can be seen from 


the fact that the residue at s = 0 represents the steady 
state of the problem. Therefore, 


y(x, 2) = V(0)t + So’ [V(r) — V(O)] X 
1 Yt, s) gat) | 
[p+ 22 Fes Sn) 7 4 
a(x, t) = fo’ [V(r) — V(0)] 
1 O(x, Sn) a) gat) | 
[2% cana” 





where Y,(a, s,) represents the partial differentiation of 
Y(a, s) with respect to s, and evaluated at s = s,, the 
nth root of the equation Y(a, s) = 0. From Eqs. (2) 
it can be seen that Y(a, s) is even in s; therefore, the 
roots s, always occur in pairs. Besides, if vibration is 
possible, all roots of the equation Y(a, s) = 0 will be 
imaginary. Let the roots be +¢éw,; then y(x, ¢) and 
6(x, t) can be put in the form 


\ 











y(x, t) = V(O)t+ fo [Vir) — V(O)] X 
Y(x, w,) COS w,(t — 7) ] 
9 
[3 1 Y,(a, Wn) wy i ” 
(6) 
A(x, t) = Sil [V(r) — V(O)} x 
[2 y O(x, wn) cos S Wy(t — t) ] in, 
- Y (a, Wy) . Wa a 





/ 


It is recognized that in Eqs. (6) Y(x, w,) and O(x, w,) are 
nothing but the mth coupled mode in natural vibration 
of the beam, with the points of attachment of the land- 
ing gear to the wing as node points. Therefore, Y(x, 
Wn), O(x, w,) and w, can be evaluated by the conven- 
tional methods used in vibration calculations (see be- 
low). 

If the quantities Y,(a, w,) can be calculated, 
y(x, t) and @(x, ¢) are completely determined. By 
differentiating Eqs. (6) with respect to ¢, the velocity 
and acceleration distributions are obtained. 


THE EVALUATION OF Y,,(a, wp) 


The differential equations (Eqs. (2)) are not easy to 
solve in a practical problem. Usually, instead of treat- 
ing a wing with continuous distribution of masses and 
moments of inertia, the masses are assumed to be segre- 
gated into several isolated points so as to reduce the de- 
gree of freedom from infinite to a finite number. In 
that case Eqs. (2) are to be replaced by a system of 
simultaneous algebraic equations of the form 

, 
V(p) = w* 2 Hoalm(a) ¥(g) + S(q)0(Q)] 
o> y 
i (7) 
O(p) = wo D dw lS(q) ¥(a) + 1(9)0(@)] 
q= 


where 7 is the number of degrees of freedom chosen in 
solving the problem. (p) or (q) denotes that the quan- 
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tity preceding that symbol is related to the pth or gth 
point; both p and g may take any value between | and 
Y. [pq is the bending influence coefficient; it gives the 
deflection at the pth point due to a unit load at the gth 
point. A,, is the corresponding torsional influence 
coefficient. Obviously, upg = Mp, ANd Age = Ngp- 

There are altogether 27 equations in Eqs. (7) with 
2r + 1 unknowns. The frequency equation is repre- 
sented symbolically by the determinant 


lax] = 0 (8) 
with 

1 
64 2 ~ ek uym(t), forl<i<r 

Ww 

l 
Cy = = wo? ca Nil (2), for r + < 1< 2r 
ay = wym(j), forl<t<r,l<j<r 
ay = \yS(j), forr +1<14< wi<j<r 
Ai = weyS(j), for 1 << i< rr+t+l<j< 2r 
ay = Agl(j), forr +1<it< 2rn,r+1<j< 2r 


When + is larger than 3, it is tedious to solve the fre- 
quency Eq. (8) directly. In that case the matrix 
iteration method is found to be handy.? By that 
method the frequency w and the mode can be found 
simultaneously. In some instances, it is easier to find 
the uncoupled bending and torsion modes first and then 
couple them together. Since the methods to solve this 
problem can be found elsewhere, no further discussion is 
given here. 

It is obvious that both the frequency and mode of 
vibration depends upon the distance 2a between the 
landing gears. Hence, Y,(x, w) can be expressed by 
the following equation: 


Y,,(x, w) = (dY/da)(da/dw) (9) 


Now, in differentiating Eqs. (7) with respect to a, one 





obtains 
dY ae [ vo) 
Hy | m(q) —— 
dO 
S@) oe < YP) + « > don 
q=1 
[m(g) ¥(g) + Sieg] f | 

40(6) yo [ rw) 

=a? ) n| SQ) =~ + 

q=1 
dO 2 “ 
1(@) 0] 2 22 000) 

(notice that \,, does not depend upon a). If regarding 


dY(p)/da, etc., as variables, by comparing Eqs. (7) and 
(10) it is found that the only condition that dY(p)/da, 
etc., are going to be finite is (for all 7 between 1 and 
2r) 


1946 


. 2 dw . du 
PP ne 2 — 
2 Aer) +04 - 

te q=l1 


[m(q) ¥(q) +S«o(@]} + L (11) 


2r 


Yo se[2ou-po 


i=r+1 ) 





where Aj, is the cofactor of a,; in the determinant 
on the left-hand side of Eq. (8). But from Egs. (7), 
if 
U(q) = m(q) ¥(q) + S(q)e() 
V(q) = Sq) ¥(q) + 1(g)0(Q) 


the following set of simultaneous equations can be de- 
rived: 


U(b) = wtm(p) HpqgU(q) + w*S(p) ae V(q) 
= q= 
(13) 


Vip) = #50) ¥ wma) + wT(P) > de VQ) 


Eqs. (13) are identical to Eqs. (7); therefore the fre- 
quency equation should have the saime roots. Now, 
the determinant that corresponds to the frequency 
equation derived from Eqs. (13) is such that its rows 
are identical to the corresponding columns in Eq. (8), 


and its columns are identical to the corresponding . 


rows. 

Therefore, for 1 < 7 € r, Ayis proportional to U(i); 
forr +1 <4 € 27, Ay is proportional to V(t — 1). 
Consequently, Eq. (11) is reduced to 


+ Up) EGrO + +e 


in 
Tue |+ 





p=1 q=l1 
. 2 dw 
> V(b) — 7 8) = 0 
p=l 
Solving for dw/da, 
d 
a > ae Ub) Ua) 
dw P= eee! =1 ¢= 
da 


” Y W@)¥@) + VP»e~) 


In the Appendix, it is shown that for symmetric load- 
ing, 


dpy,/da = R(p) + R(q) 
and for antisymmetric loading, 


dhipg _ =) ao 


fe R@ += RO) 


Therefore 
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[oe || Sxwvw | 


> | ve) Y(p) + view) | 


p=l 


= —yw> 





dw 
da 


for symmetric loading, and 


[> *ue|[ SRove | 
3 pal 





a ena 
da | sore + veer) | 


for antisymmetric loading. 
From Egs. (10), it is seen that, at x = a, 


r 


dpLag 
w? Fe lm ¥(q) + S(QO(Q)] 
a 


q=1 


adY(a) _ 
er 
In the Appendix, it is also shown that 
ditge/da = R(q) 
therefore, 
dY(a)/da = 2, RQ) U(q) 
qI= 
Consequently, the expressions for Y,,(a, w,) are found 
to be 
Y,(a, W,) = 
1 | 2, bm) 0) + 25(p) ¥(p)O() + 1(0)0%(P))] 


2 Z lm) Ve) + S(P)9)) 





(14a) 
for symmetric loading, and 


2 (a, W,) = 


dX [m(p) Y*(p) + 2S(p) ¥(p)O(p) + I(p)0?(p)] 


_ 1 ees 


. 3 =P) im(p) Y(p) + S(p)O(0)] 


p= 





(14b) 


for antisymmetric loading. In Eqs. (14), the Y’s and 
6’s are those of the mth coupled mode. 


INERTIA LOADS IN THE WING 


With Y,,(a, w,) as given by Eqs. (14), the displace- 
ment and rotation of the wing are completely deter- 
mined by Eqs. (6). Theoretically, it should be possible 
to derive the inertia loads from Eqs. (6) by differentiat- 
ing them partially with respect to x. But as Y(x, w,) 
and 6(x, w,) are known only at few isolated points, it is 
not practical to carry out the differentiation. The 
other way to determine the bending and twisting 
moments is to find the inertia loads first. Differentiat- 


ing twice the Eqs. (6) with respect to time, the linear 
and angular accelerations are obtained. 


0*y/dt? = (dV/dt) + >> Ky(t) Y(x, wa) 


(15) 
0°0/dt? = >> K,(t)O(x, wn) 


where 
2 1 SORA a 
«0 =~ lzeala -* {vo 
V(0)] cos w,(¢ — rar} (16) 


The inertia load P(p) and twisting moment 7(p) at 
station p are then 














= 0*y 070 
P(p) = m(p) = + Sip) 
0? 076 il 
a ce Meat a a 
Tp) = Sip) SF + 10) 5 


The resultant beam shears and bending moments and 
also twisting moments can be calculated by summation. 
It should be noted that, in computing the beam shears 
and bending moments, the concentrated load P(a) at 
the point of attachment of the landing gear to wing 
must be taken into consideration. P(a) is given by 
P(a) = — 2, P() (18a) 
p= 
for symmetric loading, and 


r 


P(a) = — ) 2) yp) 


p=1 


(18b) 


for antisymmetric loading. 

In Eqs. (17) it is noticed that the magnitude of 
accelerations depends upon the function K,(¢). In de- 
sign, only the maximum inertia loads are of interest. 
This requires the evaluation of the maximum of K,(¢). 
It should not present too much difficulty after V(t) is 
given. Of course, the K,(¢)’s do not attain their 
maxima at the same time. Since, in most cases, the 
maximum value of K,,(¢) decreases rapidly with increase 
in ”, it will give a reasonably conservative answer if the 
maxima of all K,,(¢) are used in Eqs. (15). 


INERTIA LOADS IN THE FUSELAGE AND TAIL SURFACES 


From the calculations in obtaining the inertia loads in 
the wing due to landing, the accelerations of the c.g. of 
the fuselage is determined. Now the fuselage can be 
regarded as an elastic beam floating in the air under 
forced vibration. The type of forced vibration is 
described by the motion of its c.g. Thus the distribu- 
tion of inertia load in the fuselage can be determined in 
a similar manner as that for the wing. The details are 
omitted in this paper. 
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In finding the inertia loads in the tail surfaces during 
landing, each surface can be regarded as a cantilever 
with its base under forced vibration. The mode of 
vibration can be obtained from the analysis of the fuse- 
lage. As a rough approximation, the mode can also be 
found directly from the analysis of the wing, by assum- 
ing that the fuselage is rigid. Because of the long arm 
from the tail to the c.g. of the airplane, very often the 
tail is subjected to high inertia loads. In some in- 
stances, these loads may cause failure if not provided 
for. 


DyNaAMIc Loap DUE To TAXIING 


The method developed in the paper can also be used 
to estimate the dynamic loads introduced from taxiing 
on a rough runway. The only change needed is the 
choice of a suitable V(t). If V(t) happens to be a 
periodic function of time and has a frequency near one 
of the coupled vibration frequencies of the airplane, the 
function K,(¢) corresponding to that mode will become 
larger and larger and finally dominate the motion. A 
state of resonance is attained, and it may become 
dangerous for the structure. 


CONCLUSIONS 


From the analysis made above, the vibrations of the 
wing excited by the impact load introduced in landing 
are found to be such that the points of attachment of 
landing gear to the wing are the node points. When the 
various modes of coupled vibration that satisfy that 
condition are available, the inertia loads can be easily 
calculated. The procedures that may be used in making 
such an estimation are listed below: 

(1) Calculate the coupled modes of vibration of the 
wing, assuming that it is supported at the landing gears. 
Usually the first three modes in each case (symmetrical 
and antisymmetrical) will be sufficient in estimating 
the inertia loads. 

(2) Calculate the maximum values of K,(¢) from Eq. 





1946 


(3) Calculate the accelerations at all stations by 
Eqs. (15), using the maximum values of K,(¢) obtained 
from step (2). 

(4) Calculate the inertia loads at the stations by 

(17) and the reaction of the landing gear by 


Eq. (18). Sum up these inertia loads to find the 
beam shears, bending ‘moments, and twisting mo- 
ments. 


A large part of the labor will be spent in finding the 
modes and frequencies of the coupled vibration. Be- 
cause of the fact that these vibrations have their node 
points at the points of attachment of the landing gear to 
the wing, they, in general, are quite different from those 
of afree wing. Therefore, that part of the labor seems 
to be unavoidable if an accurate answer is required. In 
summing up the inertia loads and moments, it is advis- 
able to calculate those due to each mode separately and 
then add them together. This is necessary because the 
maximum value of K,,(¢) is used in step (3). 

The accuracy of this method, of course, depends upon 
the accuracy of the functions V(a, t) and V(—a, ?). 
These functions, as stated in the paper, can be estimated 
from the experimental data of drop test of the landing 
gear in question. In the conventional drop test, a 
rigid mass is attached at the top of the landing gear. 
The result thus obtained can be regarded only as an 
approximation. If there is a better way to figure out 
the magnitude and time history of these functions, the 
method outlined in the paper should give a close 
estimation of the dynamic loads introduced during 
landing. Therefore, it is highly desirable to have 

some work done in that respect. It will also be helpful 
if experimental data are available to check the accuracy 
of this method. 
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Appendix 


INFLUENCE COEFFICIENTS IN BENDING AND TORSION 


’ The influence coefficients in bending can be easily 
obtained by the conjugate beam method, The results 
are summarized in the following: 


(1), Symmetric Bending (Fig. 2a) 


(i) p<aq<a 
(a — x)(x — p)dx 


— — p)R wr 
Mp (a — p)R(q) J B 





(ii) p<a,qg>a; orp>a,q<a 





e=(a- p)R(q) 
(iii) p>a,qg>a 


Mpg = (a — p)R(g) + fi (a — 2e — p)dx _ 


-of (@« - =f 
where R(q) is given by 


R@ =(@-9 [8+ [e=8 











Fic. 


for. 


(iii) 


Mpg ’ 


whe 


for « 


for ¢ 








ns by 
tained 


ms by 
sar by 
id the 
> mo- 


ng the 

Be- 
r node 
rear to 
| those 
seems 
1. In 
advis- 
ly and 
se the 


| upon 
‘a, t), 
nated 
nding 
at. a 
gear. 
aS an 
e out 
3, the 
close 
uring 
have 
‘pful 
racy 


Engi- 
York, 


ds in 
Néw 


DYNAMIC LOADS IN STRUCTURES DURING LANDING 387 





(49152 Womb sors 
Pie = 35 
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Fic. 2. Loading conditions for finding bending influence coeffi- 


cients. (a) Symmetrical; (b) antisymmetrical. 


for gq < a, and 
R(q) = (a — q) fo* (dx/B) 
forg > a. 


(2) Antisymmetric Bending (Fig. 2b) 


4 f° Ga se = pie 
a p B 


orp >a,q<a 
= (a — p)R@) 


fi) p<a,q<a 





yg = (a — p)R(Qq) 


(ji) p>a,q>a 
Gh ptig + fem aie 


oo f e=pe 
where R(q) is given by 


a— * x*dx * x(a — 
pe a ee 
for gq < a, and 
R(q) = 
for g < a. 








[(a — g)/a*] f° (x*dx/B) 


The geometric meaning of R(q) is that it is the slope 
at the support due to unit loads at g. The magnitude 
of (du/da)pg can be found by direct differentiation. In 
symmetric bending it is found that 


duy»,/da = R(p) + RQ) 
and in antisymmetric bending, 
dpy,/da = (g/a)R(p) + (p/a)R 


It is noticed that the bending coefficients listed above 
can be expressed in terms of three fundamental in- 
tegrals—namely, f;° dx/B, Js (a — x)dx/B, and 
JS;°(a — x)*dx/B. These integrals can be evaluated 
numerically. In carrying out the numerical integra- 
tion, the intervals used should be smaller than the 
distance between the stations chosen in the vibration 
analysis. If r stations are used in the analysis, 2r 
intervals should be used to evaluate these integrals. 
The whole work can be done in one suitable table. 

As can be seen in the body of the paper, the supports 
of the beam have nothing to do with the twisting of the 
beam. Therefore, the influence coefficients in torsion 
are the same as those for a free wing. The formulas are 
collected below: 


(i) p< 4, 

Ang = to + So? dx/C 
(ii) p>q, 

Noe = do + Jot dx/C 
where for symmetric torsion, 


do = — LI (Apy — %0)/DI 


r 


and for antisymmetric torsion, 


ho = 0 
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